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^ _  _ - _ SUMMARY  •/ 

0  i  to  t’dv  ■'  d  j 

The  Naval  Ocaan  Systems  Center  (NOSC) ,  in  collaboration  with  tha  Canadian 
Dafanca  Raaaarch  Batablishaant  j  Pacific  (DRIP),  conducted^*  Joint  Ocaan  Nava 
investigation  Pro j act  ( JOMIP) { to  avaluata  tha  datactability  of  ocaan  wava 
structures  on  imaging  synthetic  aperture  radar  (SAR) .  This  pro jact  <made?  use  l -> 
of  relatively  wall  understood  ^Kelvin  surface  ship  wake  patterns  generated 
under  controlled  and  wall  documented  surface  environmental  conditions  to 
isolate  SAR  image  parameters. 

SEASAT  SAR  images  of  the  ocaan  surface  established  the  importance  of 
wakes  in  ship  surveillance  from  space.  -Th^  jjse  of  waves  of  known  wavelength 
and  direction  provide  the  opportunity  to  evaluate  the  SAR  contribution  to 
ocean  wave  forecasting. \  Airborne  SAR  operated  by  the  Canadian  Centre  for 
Remote  Sensing  ^CCRS )  gathered  the,  doppler  histories  and  high  resolution 
images  were  generated  by  the  Environmental  Raaaarch  Institute^  of  Michigan . 


<1— SEASAT- like  images  made  with  L— band  SAR  are  presented:  ( 1 )  to  suggest  how 
ship-generated  surface  roughness  combines  with  the  velocity-bunching  mechanism 
in  calm  ocean  areas  to  produce  the  unusually  narrow  wakes  observed  for 
azimuth- traveling  ships  on  SEASAT  images,'  and  (2)  to  estimate  conditions  under 
which  SAR  image  modulation  mechanisms  can  be  expected  to  produce  wake  Images. 

A  method  is  described  for  using  the  Kelvin  transverse  ship  wake  wave 
component  to  quantitatively  evaluate  the  contributions  of  various  SAR  ocean 
wave  Imaging  mechanisms.  method^ makes  use  of  the  narrow  sector  of  sur¬ 

face  roughness  generated  by  a  ship  along  its  track  to  produce  SAR  image’s  of 
the  longest  waves  in  its  wake  system  on  flat  calm  water,  k  By  varying  ship 
headings  and  speeds,  relative  SAR  image  modulation  mechan: sms  can  be  quan¬ 
titatively  evaluated  as  a  function  of  surface  wave  propagation  direction, 
wavelength,  radar  incidence  angle,  and  to  some  extent,  height.  under  near 
aero  local  wind  conditions,  ship-induced  roughness  may  also  be  used  to 
determine  ocean  swell  conditions.  f 


BACKGROUND 


The  existence  of  ship  wakes  on  SEASAT  SAR  images  has  been  discussed  by 
multiple  authors  { Fn  and  Holt,  1982;  Vesecky  and  Stewart,  1982),  and  the  use 
of  Kelvin  surface  wake  components  for  investigating  surface  wave  SAR  image 
modulation  mechanisms  has  been  pursued  by  NOSC  for  several  years  (Hammond,  et 
al.,  1983).  Che  area  of  specific  interest  has  been  the  narrow  wake  angles 
observed  on  many  SEASAT  images .  Suggested  responsible  mechanisms  for  these 
narrow  wake  images  include  (1)  increased  radar  roughness  associated  with  ship¬ 
generated  surface  wavelengths  near  the  Bragg  reinforcement  length  («30  cm  for 
SEASAT);  (2)  enhanced  slopes  associated  with  these  shorter  wavelengths;  (3) 
the  surface  manifestations  of  ship-generated  internal  waves;  or  (4)  some 
combination  of  the  above.  Controversy  over  wakes  is  inevitable  when  we  try  to 
explain  specific  wave  observations  with  imaging  theory  derived  with  statis¬ 
tical  ensemble  characteristics  as  variables.  The  application  of  SAR  to  naval 
ocean  surveillance  missions  reguires  detailed  understanding  of  both  discrete 
wake  patterns  and  statistical  wave  ensembles. 


A  SEASAT  SAR  data  base  including  49  ship  wakes  over  the  full  range  of 
ship  headings  was  compiled  from  JPL  images  which  were  processed  digitally  and 
optically.  This  data  base  produced  the  statistics  shown  in  table  1  as  a 
function  of  ship  track  orientation  relative  to  the  SAR  track.  It  should  be 
pointed  out  that  an  estimate  of  ships  heading  was  prereguisite  for  inclusion 
in  this  data  base.  Because  the  resolution  of  these  images  varied  from  25  to 
40  meters  and  the  longest  wave  in  a  "typical"  (Webster,  1982)  15-knot  ship 
wake  is  only  38  meters  the  expected  wake  image  included  only  the  common  Kelvin 
wake  envelope  arms  at  approximately  VI  9.5  degrees  and  perhaps  an  elongated 
turbulent  patch  extending  for  some  distance  directly  aft  of  the  ship. 

Tfcble  1 .  SEASAT  ship  wake  image  data  summary.  . 
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While  each  of  the  49  ships  could  display  some  combination  of  these  three 
features,  table  1  indicates  that  35  percent  of  the  ships  displayed  neither 
port  nor  starboard  envelope  arms.  The  turbulent  wake,  usually  an  absence  of 
radar  return  against  a  saturated  background,  was  the  most  prevalent  and  sur¬ 
prisingly  persistent  wake  feature  averaging  8  km  long.  The  turbulent  wake  was 
the  only  feature  observed  for  ships  traveling  in  the  radar  range  direction. 
As  ships  headings  became  more  along  track,  wake  arms  began  to  appear,  usually 
with  the  radar  range  traveling  wake  component  side  of  the  wake  first,  and  both 
sides  became  increasingly  common  as  the  relative  track  orientation  angle 
approached  zero.  Wake  arms  averaged  5  km  long  on  SEASAT  images  but  their 
included  angles  varied  widely  from  the  expected  4^1 9.5-degree  angle  associated 
with  the  Kelvin  envelope  where  the  highest  amplitude  waves  are  known  to  be 
located  (mui,  1962). 

The  distribution  of  observed  wake  half-angles  (each  side  considered 
separately)  displayed  the  two  distinct  modes  shown  in  figure  1.  The  upper 
mode  peaks  near  the  anticipated  Kelvin  envelope  angle,  but  a  second  unexpected 
peak  near  3  degrees  is  also  a  dominant  feature. 


Figure  1.  Observed  distribution  of  SEASAT  ship  wake  angles. 


During  high  surface  winds  with  highly  nonlinear  imaging  conditions 
(Alpers,  1983),  no  wake  features  were  observed.  With  lower  winds  the 
turbulent  wake  appears  first  as  a  reduced  radar  return  against  a  saturated 
background.  As  wind  speed  approaches  the  linear  imaging  regime  («  3.6  m/s  for 
SEASAT  with  four  resolution  cells  per  ocean  wavelength  [Vhsecky,  et  al., 
1982]),  smeared  lines  may  appear  on  high  resolution  images  along  the  Kelvin 
envelope  arms  (range-traveling  arm  first)  and  both  arms  may  become  apparent 
for  range-traveling  ships  where  the  wake  component  waves  propagate  in  the 
linear  8AR  modulation  domain  (Alpers  et  al.,  1981).  At  still  lower  wind 
speeds  there  is  insufficient  radar  cross-section  <o Q>  for  the  wake  component 
waves  to  modulate  the  image  and  these  waves  gradually  disappear  (upper  mode  in 
figure  1)  and  are  replaced  by  the  narrow  wake  (lower  mode).  We  hypothesized 
these  smaller  observed  wake  angles  were 


(1)  dependent  on  short  gravity  wave  surface  roughness  generated  by  the 
ship,  and 

(2)  observable  only  under  relatively  calm  local  wind  conditions. 

During  the  summer  of  1983,  the  Naval  Ocean  Systems  Oenter  (NOSC)  in 
conjunction  with  the  Canadian  Defence  Research  Establishment,  Pacific  (DREP) 
embarked  on  a  joint  project  to  observe  ship  wake  structures. 

the  U.S.  Navy's  objectives  in  the  U. S. /Canadian  Joint  Ocean  Wave 
Investigation  Project  (JOWIP)  were  to 

(1)  resolve  the  narrow  wake  image  controversy 

(2)  test  ocean  surface  wave  imaging  theories  using  the  well  known 
components  of  the  Kelvin  wake  pattern 

(3)  evaluate  SAR  detection  of  ship-generated  internal  waves  (IW) 
relative  to  naturally  occurring  IW  patterns. 

DREP,  which  is  a  highly  capable  Canadian  Forces  laboratory  in  Esquimalt, 
British  oolunbia,  shares  similar  interests  with  NOSC.  Because  of  this  and 
DREP's  considerable  experience  with  ship-generated  internal  waves  and  access 
to  the  Canadian  centre  for  Remote  Sensing  (CCRS)  airborne  simultaneous  dual 
polarization  X-L  band  SAR,  a  joint  venture  sharing  the  surface  truth  collec¬ 
tion  burden,  instruments,  and  data  seemed  appropriate  to  SAR  imaging  problems. 

One  serious  problem  in  any  type  of  SAR  evaluation  has  been  the  paucity  of 
image  data  obtainable  for  the  multitude  of  uncontrollable  environmental, 
flight,  and  radar  variables  involved.  Our  plan  was  to  concentrate  on  the 
SEASAT  parameters  (L-band,  horizontally  polarized, ~23°  surface  incidence 
angle)  and  obtain  high  confidence  data  on  a  restricted  set  of  variables  before 
extending  the  cases  to  be  considered. 

Since  statistical  description  of  the  ocean  surface  waves  injects  large 
uncertainties,  we  used  the  relatively  well  known  wavelengths  and  propagation 
directions  within  the  Kelvin  surface  wake  of  a  large  dedicated  U.S.  Navy  ocean 
tug.  Ihe  plan  called  for  demonstrating  reproducible  results  under  conditions 
critical  for  verification  of  current  image  modulation  theory;  first,  on  flat 
calm  waters  where  the  scattering  waves  would  be  generated  by  the  ship  and 
second,  on  a  modestly  wind- roughened  surface  where  background  waves  would  pose 
minimum  interference  in  resolving  the  longer  ship  wake  components.  The  ship 
was  to  travel  (1  )  in  the  same  direction  as  the  radar  orienting  its  shorter 
wake  components  in  the  radar's  range  direction,  the  situation  which  produced 
narrow  wake  angles  on  SEASAT  images;  and  (2)  perpendicular  to  the  radar  track 
to  orient  the  longer  wake  wavelengths  in  the  range  direction. 

In  order  to  separate  surface  and  internal  wave  effects  it  became  neces¬ 
sary  to  operate  in  two  separate  areas  sheltered  from  wind  and  swell  as  much  as 
possible  —  a  primary  location  with  as  little  stratification  as  possible,  and 
the  other  location  highly  stratified.  Choice  of  actual  areas  involved 
tradeoffs  against  ship  transit  times,  Canadian  aircraft  flight  times,  and 
basing  requirements  as  well  as  local  wind,  swell,  and  internal  wave 


conditions.  After  all  factors  were  considered,  we  selected  our  primary  sites 
in  the  HUMES  range  in  Dabob  Bay,  Washington  and  a  very  highly  stratified 
location  in  Knight  inlet  in  British  ODlunbia.  Surface  ship  wake  tests  were 
conducted  in  Dabob  Bay  between  26-28  July,  and  the  IW  generation  tests  were 
conducted  in  Knight  Inlet  on  2-3  August  1983* 

The  surface  wave  experiments  involved  higher  speed  (12-16  knots)  ship 
runs  on  Dabob  Bay.  Results  of  these  tests  are  reported  here  by  HOSC  and  the 
results  of  the  ship-generated  and  natural  internal  wave  observations  will  be 
prepared  by  DREP.  Ship-generated  internal  waves  were  produced  by  both  the  USS 
pUAPAW  and  CFAV  ENDEAVOUR  running  at  1  -5  knots  in  the  highly  stratified  waters 
of  Knight  inlet,  The  SAR  images  of  these  wakes  revealed  weak,  dark  (absence 
of  return)  lines  against  a  wind-roughened  surface.  Ho  corresponding 
subsurface  evidence  of  internal  wave  wake  activity  was  detected  on  Dabob  Bay 
during  repeated  CFAV  ENDEAVOUR  echo so under  passes  through  the  QUAPAW  wake 
system  (Hughes,  1983).  As  will  be  seen  later,  the  absence  of  wind  roughening 
during  the  Dabob  SAR  runs,  which  resulted  in  narrow  wake  images,  removed  the 
SAR  clutter  background  against  which  internal  waves  are  normdly  seen. 

The  Canadian  Centre  for  Remote  Sensing  aircraft  SAR  flew  1 5  passes  over 
the  Dabob  location  with  the  US  Navy  ocean  tug  (USS  QUAPAW)  laying  down  wake 
patterns.  The  Canadian  oceanography  ship,  CFAV  ENDEAVOUR ,  collected  wind  and 
surface  wave  height/slope  data.  A  stereo-photo  fixed-wing  aircraft  and  a 
helicopter  recorded  the  events  on  vertical  and  oblique  photographs,  respec¬ 
tively.  The  SAR  imaging  geometry  was  held  as  close  as  possible  to  that  of 
SEASAT,  with  primary  emphasis  on  L-band  HH  polarization,  and  the  image  times 
were  selected  to  the  maximum  extent  possible  for  wind  and  wave  conditions. 
The  CCRS/ERIM  SAR  produces  slant  range  and  azimuth  resolutions  of  3  and  5 
meters,  respectively,  for  stationary  targets  using  a  coherent  integration  time 
of  1.4  seconds.  Even  after  translation  to  the  ground  plane  and  correction  for 
motion-induced  image  smear,  a  16-meter  maximum  effective  resolution  was 
achieved  when  estimated  by  the  equations  of  Alpers  et  al.,  (1981).  It  was 
hoped  that  this  SAR  would  resolve  the  longer  transverse  waves  (  «  38  m)  in  a 
15-knot  ship  wake  under  proper  surface  conditions  using  the  effective 
resolution  cells  per  ocean  wavelength  (n  _>  2)  criterion  of  Vesecky,  et  al; 
1982.  * 


SHIP  HAKES  AND  SAR  IMAGE  PARAMETERS 


The  following  brief  background  on  ship  wake  structures  and  SAR  imaging 
geometry  will  provide  the  basic  equations  and  parameters  needed  for  the 
analysis  which  follows . 

SURFACE  HAKE  STRUCTURES 

A  ship's  hull  imparts  velocity  to  the  water  at  all  angles  between  0  and 
90  degrees  to  the  ship  track.  The  initial  amplitudes  depend  on  the  shape  of 
the  hull  but  outside  the  immediate  disturbance  area  the  resulting  "free  waves" 
of  wavelength  (1)  and  velocity  (c),  propagating  in  any  direction  (<*)  are 
related  to  the  ship's  speed  (V)  (Havelock,  1934)  by: 

a  -  cos  *  C/V  ( 1 ) 

Ha)  -  | -  V2  cos2  a  (2) 

He  may  summarize  the  velocity  vectors  originating  from  a  ship  at  point  P  as 
shown  in  figure  2. 


Figure  2.  Hake  velocity  vectors  as  a  function  of  propagating  angle  (a). 


The  energy  of  the  waves  in  direction  a  originating  at  point  P  with  phase 


velocity  C  will  propagate  with  group  velocity  C  *  C  /2  as  indicated  by  the 
P  9  P 


smaller  semicircle  in  figure  2.  Note  that  although  the  deep  water  wave 

velocity  (Cp  -  >/gX/2w )  will  depend  on  ship  speed  (V),  the  ratio  C^/V  for  any 

specific  a  is  constant.  This  constant  expansion  rate  agrees  with  our  obser¬ 
vation  that  ship  wake  structures  expanding  in  the  cross  track  direction  appear 
stationary  with  respect  to  the  ship  in  photographic  and  SAR  images. 


Since  the  common  measurable  parameters  on  ship  wake  images  are  ( 1 )  the 
ship  heading  (H^),  (2)  the  vee-angle  or  commonly  half  vee-angle  (0  ) ,  and  (3) 
one  or  two  wake  envelope  arm  lengths,  it  is  advantageous  to  express  wake 


component  parameters  for  all  points  (Q)  within  the  envelope  arms  in  terms  of 
the  polar  coordinates  (0  ,r)  with  origin  approximately  at  the  ship  as  shown  in 
figure  3. 


Figure  3.  Hake  envelope  coordinates  (X,  0). 


This  geometry  allow*  the  solution  for  wake  coordinate  8  in  terms  of  propaga¬ 
tion  directions  (a)  independent  of  ship  speed  and  hull  shape,  as  given  by 
equation  (3)  which  associates  two  propagation  directions  (a)  with  each 
observed  wake  coordinate  angle  ( 8Q ) . 


a 


arc tan 


t  t  ( 1  -  8  tan2  6 ) 
4  tan  8 


(3) 


The  longer  wavelength  Kelvin  wave  components,  which  move  more  or  less  in  the 
direction  of  ship,  are  termed  transverse  (aT  <  35*)  and  those  more  or  less 
crosstrack  propagating  waves  (35*  <  oD  <  90*)  are  termed  divergent.  Havelock 
(1934),  Inui  (1962),  and  Lighthill  (1978)  in  their  pioneering  work  on  the  wave 
resistance  of  ships,  have  shown  that  it  is  possible  to  obtain  a  good  approxi¬ 
mation  to  surface  wake  amplitude  A(X,  8)  for  any  point  in  the  wake  in  terms  of 
a  linear  combination  of  the  amplitudes  ^(o^)  and  corresponding  to 

these  two  "elementary  waves."  This  method  is  applicable  to  the  entire  range 
of  ocean  vessels,  Froude  numbers  -  V/^/gL  <  0.35  (Inui,  1962),  for  locations 
further  than  several  ship  lengths  aft.  Variations  of  these  equations  can 
accommodate  all  points  within  the  outer  edge  of  the  Kelvin  vee-wake  pattern 
(8  ■  19.5*,  a  -  35)  where  the  familiar  larger  amplitude  visible  "cusp"  waves 
mark  the  transition  between  transverse  and  divergent  wake  components. 


Assuming  that  the  SAR  wake  pattern  on  a  flat  uniformly  Bragg  roughened 
aur face  results  from  surface  wake  disturbances,  the  wavelengths,  propagation 
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directions,  amplitudes,  and  resulting  slopes  can  be  obtained  for  resolution 


elements  along  any  observed  wake  arm  (6  -  0Q).  Propagation  directions  (a)  for 

a  common  SAR  image  wake  coordinate  (0  ■  3.5s)  are  a  «  3.5*  and  a  *  83*. 

0  T  D 

The  corresponding  wavelengths  for  the  transverse  (a^  -  3.5*)  and  divergent  (aQ 
»  83*)  components  are  shown  in  table  2  for  various  ship  speeds. 


Table  2.  Surface  ship  wake  component  wavelengths  along  the  0Q  «  3.5* 
wake  coordinate  angle  for  various  ship  speeds 


Ship  Speed  (knots) 

12. 

I! 

12 

H 

15 

Divergent  Wavelength  (m) 

0.25 

0.37 

0.43 

0.50 

0.57 

Transverse  Wavelength  (m) 

16.9 

24.4 

28.6 

32.2 

38.1 

JOWIP  ship-wake-component 

amplitudes 

several 

ship 

lengths  aft 

were  i 

mated  to  be  between  5  and  25  centimeters  before  the  experiment.  For  planning 
purposes,  we  used  10  centimeters  1000  meters  aft.  Preliminary  stereo  photo¬ 
graphic  measurements  and  shipboard  wave  height  time  series  data  have  confirmed 
this  but  precise  variation  with  propagation  direction  remains  to  be  extracted 
from  wake  photography.  The  length  of  time  that  short  gravity  divergent  waves 
(  ~30  cm)  could  survive  in  the  open  ocean  is  another  area  of  uncertainty.  In 
calm  ocean  areas,  where  the  population  of  wind-generated  short  wavelengths  is 
depleted  and  strong  interactions  minimized,  something  like  the  10-minute 
energy  e-folding  time  estimated  for  viscous  damping  (Lighthill,  1978)  might  be 
appropriate.  Even  this  short  interval  could  yield  a  4.5-kilometer  wake  for  a 
15-knot  ship. 

INTERNAL  WAVES 


Since  the  se  (ration  of  Kelvin  wake  and  internal  wave  effects  on  SAR 
images  is  the  prim.  'Mective  of  the  NOSC  part  of  the  JOWIP  experiment,  it 
was  desirable  to  choci  measurement  site  where  surface  ships  were  not  able 
to  produce  internal  wav..  Depending  principally  on  ship  speed  and  density 
structure  of  the  water  cox^n,  two  diBtinct  physical  mechanisms  associated 
with  surface  ship  motion  are  capable  of  generating  internal  wave  wakes:  water 
vertical  displacement  by  the  hull  and  collapse  of  the  turbulent  wake. 

Hull-induced  internal  waves  having  wavelengths  near  the  ship  length  can 
be  significant  if  the  ship  speed  is  matched  to  their  phase  velocity  in  the 
medium.  This  phase  velocity  is  given  by  a  dispersion  relationship  which 
depends  on  the  density  gradient  in  the  water  column.  These  waves  also  exhibit 
a  maximum  frequency  which  is  given  by  the  quantity  known  as  the  Brunt-Vaisala 
frequency  ( N { z ) ) . 


N(z) 


q  dp(z)  j  1//2 
p(z)  dz  J 


radians  sec 


-1 


(4) 


SAR  TRACK 


SAR  WAKE  OBSERVATIONS 


Figures  10,  11,  and  12  show  representative  comparisons  between  SEAS AT  and 
the  CCRS/ERIM  SAR  wake  images  for  azimuth-,  oblique-,  and  range- traveling 
ships.  Fbr  the  JCMIP  L-band  SAR,  strong  narrow  wake  modulations  were  observed 
as  expected  for  azimuth-traveling  ships  under  low  wind  conditions  (~2.0  m/s). 
Cosparison  with  SEASAT  generally  showed  similar  backgrounds,  contrasts,  and 
wake  angles.  Inner  wake  detail  is  usually  absent  from  SEASAT  images  because 
wake  wavelengths  ( ~40  m)  are  short  when  compared  to  the  processed  resolutions 
(  ~25  m).  Wake  angles  for  most  of  the  CCRS/ERIM  SAR  images  appear  narrower 
because  these  optically  correlated  images  are  in  the  slant  rather  than  the 
ground  range  plane . 

In  order  to  resolve  ship  wake  angles  on  SAR  images,  the  wake  arm  must 
exceed  the  theoretical  miniraun  lengths  shown  in  table  7.  The  lengths  are 
based  on  only  one  resolution  cell  between  the  wake  centerline  and  the  wake 
arm,  ideal  SAR  performance,  and  no  error  in  measurement.  Note  how  the  higher 
JOWIP  range  resolution  reduces  the  required  wake  lengths.  From  the  first  two 
entries  in  table  7  it  is  apparent  that  wake  angle  measurement  errors  of  0.5  to 
1  degree  are  to  be  expected  (observed  wake  arm  lengths  roughly  2  and  7  km  in 
figure  10)  even  under  ideal  conditions. 

Table  7.  Minimum  wake  arm  lengths  for  SAR 
wake  angle  measurements  (23®  incidence). 


Wake  Half  Angle  0Q  (deg) 

JOWIP  (km) 

SEASAT  (km) 

CV  990  (km) 

0.5 

1.84 

5.75 

4.83 

1.0 

0.92 

2.88 

2.42 

1.5 

0.61 

1.90 

1.60 

2.0 

0.46 

1.43 

1.20 

2.5 

0.37 

1.15 

0.97 

3.0 

0.31 

0.95 

0.80 

3.5 

0.27 

0.83 

0.70 

4.0 

0.23 

0.73 

0.61 

In  the  upper  SEASAT  image  in  figure  10A  for  azimuth  traveling  ships,  wake 
angles  for  the  bottom  center  wake  measured  between  3  and  4  degrees.  In  the 
slightly  wind-roughened  area  half  way  up  the  image,  one  wake  arm  (probably 
upwind)  is  missing.  7he  lower  CCRS/ERIM  image  (figure  10B>  shows  a  combina¬ 
tion  sequence  where  the  positive  return  wake  arms  disappear  almost  immediately 
and  are  replaced  by  a  dark  (absence  of  return)  linear  scar  where  the  ship 
passed  through  the  locally  wind- roughened  water  surface  area.  This  wind 
roughening  was  simultaneously  recorded  photographically  and  is  estimated  to 
correspond  to  wave  amplitudes  of  a  few  centimeters.  The  wake  modulation  is 


Table  6.  Da bob  Bay  surface  wind  suaaary. 


SAR 

Fit -Pa a 
(nuaber) 

Tlae 

Interval 

Speed 

(a/a) 

Direction 
( °T) 

Tine 

Interval 

(PDT) 

Speed 

(a/a) 

Direction 

CT) 

1-4 

No  data 
available 

5-3 

1024-1025 

5.1  ♦  0.3 

186-220 

1025-1031 

5.8  +  0.5 

193-219 

5-4 

1103-1108 

7.7 

205 

6-1 

1518-1523 

6.0 

210 

6-4 

1643-1648 

6.4  +  0.9 

180-199 

1652-1658 

5.3  +  0.5 

180-219 

7-1 

1605-1610 

11.2 

235 

7-3 

1706-1711 

9.7 

230 

7-5 

1822-1826 

7.0 

235 

7-7 

1932-1937 

6.5 

235 

8-1 

0724-0727 

2.2  +  0.2 

239-253 

0735-0739:40 

1.8  ♦  0.1 

214-240 

8-2 

0812-0817 

1.5  +_  0.3 

180-230 

0819-0826:30 

2.0  +  0.1 

180-216 

8-3 

0859-0904 

2.5  _+  0.3 

189-206 

0906-0910:20 

2.8  ♦  0.2 

188-204 

8-4 

0943-0948 

2.5  +_  0.3 

191-226 

0944-0955:30 

2.0  +  0.2 

198-224 

9-1 

1313-1317 

<1 

SW 

1317-1324 

<1 

SW 

9-2  1353:30-1354: 

30  <1 

SW 

1356-1401:10 

<1 

SW 

9-3 

1442-1447 

<2 

SW 

1449-1445:15 

<2 

SW 

NOTES: 

Referenced  to  10  aetera  above  surface 

Bi  tries  where  the  ENDEAVOUR  wake  sampling  interval  is  emitted  are  rough  esti¬ 
mates  since  digital  records  have  not  been  processed. 


Threshold  probleas  for  very  low  speed 


SAR  AMD  SURFACE  TRUTH  DATA  OBTAINED 


During  the  3  days  of  the  surface  ship  wake  phase  of  the  experiment,  SAR 
coverage  of  Dabob  Bay  produced  15  SAR  images  during  the  time  intervals  and 
with  the  ship  speeds  and  headings  swmarized  in  table  5.  The  intended  ship 
speed  was  12  knots  through  run  8-1  when  it  was  increased  to  15  knots.  The 
ship  speed  entries  in  table  5  are  presented  to  one  decimal  place  where  they 
are  estimated  and  to  two  decimal  places  in  the  cases  where  they  are  accurately 
known. 


The  SAR  image  times  were  coordinated  to  fall  within  the  interval  of  the 
QUAPAW  high  speed  run  as  shown.  Most  of  the  runs  were  made  with  both  QUAPAW 
and  the  SAR  aircraft  on  south  to  north  tracks.  Notable  exceptions  to  this 
were  passes  7-1,  7-3,  and  9-1  with  perpendicular  tracks  and  pass  9-2  where  the 
platforms  ran  in  opposite  directions.  The  SAR  image  from  pass  1-4  was  col¬ 
lected  during  a  SAR  pass  over  the  SAR  calibration  area  and  surface  truth  is 
not  available. 

An  attempt  was  made  to  compile  data  for  each  selected  SAR  image  for  a  5- 
minute  background  interval  as  well  as  during  the  time  the  DREP-instrvmented 
CFAV  ENDEAVOUR  traversed  the  wake  pattern.  The  surface  data  taken  by  DREP 
from  the  CFAV  ENDEAVOUR  includes  extensive  records  of  wind,  wave  height,  wake 
slope,  echo-sounder,  CTD  profiles,  and  precise  (approx  _+  2  meter)  latitude- 
longitude  tracks  of  pUAPAW,  ENDEAVOUR  and  Cessna  vertical  photo  aircraft.  The 
bulk  of  this  time  series  data  is  available  on  VAX/VMS  1600  BPI  9-track 
computer  compatible  digital  tapes  prepared  by  DREP  and  is  available  through 
NOSC.  These  data  have  been  carefully  corrected  to  geocoordinates,  calibrated, 
and  extensively  cross  checked.  Appendix  A  summarizes  the  available  digital 
data. 


Appendix  B  summarizes  the  times  that  vertical  photographs  were  taken 
within  the  wake  system.  Wake  coordinates  associated  with  each  of  these 
exposures  are  being  determined  from  the  DREP  digital  positioning  data  before 
selecting  specific  frames  for  analysis.  Resolution  from  synchronized  stereo 
pairs  were  used  to  take  vertical  photographs  throughout  the  wake  system  to 
measure  the  lengths  and  heights  of  the  anticipated  short  gravity  waves.  The 
fixed  baseline  stereo  images  are  on  70  mm  film,  nominally  60  meters  square 
with  75  percent  overlap  and  2  cm  ground  resolution.  A  selected  example  (run 
8-3)  where  rough  track  reconstruction  could  be  made  from  known  objects 
included  in  the  images  will  be  UBed  in  the  analysis  presented  here. 

Appendix  C  provides  the  times  associated  with  70  mm  oblique  photos  taken 
by  a  roving  helicopter  from  approximately  100  feet  throughout  the  experiment. 
The  photographs,  which  provide  some  of  the  best  qualitative  estimates  of  the 
overall  surface  roughness  patterns  in  the  SAR  images ,  will  also  be  used  in  the 
data  analysis. 

Because  of  its  critical  importance  to  water  surface  radar  scattering, 
background  winds  are  presented  in  table  6  for  the  time  intervals  and  SAR 
images  shown. 

With  this  background  information  it  is  now  possible  to  suamarize  per¬ 
tinent  observations  in  the  SAR  images  obtained. 
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Figure  8.  Velocity  bunching. 

(a)  Surface  elevation  associated 
with  ocean  wave  traveling  in 
azimuthal  direction,  (b)  Azimuthal 
shift  of  scatters  in  image  plane, 
x  is  azimuthal  axis  on  real  surface 
and  x  on  mapped  surface,  i.e,  in 
image  plane.  Shift  x  -  xq  is  given 
by  velocity  of  scatterers  in  range 
direction,  (c)  Relative  power 
(intensity)  per  azimuthal  resolution 
cell  as  function  of  azimuthal 
position  x  for  different  values  of 
c  defined  By  equation  (11)  (after 
Alpers  and  Rufenach,  1979) 


Figure  9.  The  dependence  of  non- 
dimensional  modulation  transfer 
functions  describing  tilt,  weak 
hydrodynamic,  and  velocity  bunching 
modulation  on  azimuthal  angle  (after 
Alpers  et  al.,  1981). 


Table  4.  Estimated  maximum  crest- to- trough  image  modulation. 


SAR  Platform 

SEASAT  (X  - 

100  m,  A  ■  0.5  m) 

JOWIP  (X  - 

35  m,  A  ■  0.1  m) 

Ship  Travel 

Range 

Azimuth 

Range 

Azimuth 

Tilt 

2.4  dB 

0.2  dB 

1.4  dB 

0.1  dB 

Hydrodynamic 

1.1 

0.0 

0.7 

0.0 

Velocity  Bunching 

o 

e 

o 

11.0 

0.0 

7.1 
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ROC)  -  >4.5  K  u 


,  2  . 
•in  4, 


(10) 


L  -  +  K  J 

where  u"1  is  termed  the  relaxation  time  constant.  Following  standard  practice, 
we  let  p  “  0  and  estimate  the  maximum  RCS  modulation  contribution  as  Ac  “  <0g> 
[1  _+  4.5  K  A]  sin  4  (Vesecky  and  Stewart,  1982).  For  the  above  cases  of 
100-meter  and  35>meter  range- traveling  waves,  this  gives  a  maximal  crest- to- 
trough  modulation  of  1.1  dB  and  0.7  dB,  respectively.  The  first  order  esti¬ 
mated  hydrodynamic  modulation  is  zero  for  all  azimuth- traveling  waves  (•  ■  0). 
Since  magnitudes  given  by  equation  (10)  are  thought  to  be  low  (Wright  et  al., 
1980),  additional  effects  such  as  nonlinear  wave  interactions  or  other 
scattering  mechanisms  need  to  be  investigated. 

VELOCITY  BUNCHING 

When  long  wave  orbital  motions  have  spatially  varying  radial  velocities, 
the  doppler  coordinate  is  misaddressed  causing  nonuniform  asimuth  (X)  dis¬ 
placement  of  scattering  elements  in  the  image  plane.  This  occurs  for 
uniformly  roughened  surfaces  independently  of  other  RCS  modulation  mechanisms 
and  can  accentuate  or  even  introduce  fictitious  wavelengths.  To  linearly  map 
wave  fields  into  SAR  images  and  preserve  ocean  wave  spectra  requires  (Alpers 
and  Rufenach,  1979)  that, 

c  -  1*^1  A  <  0.3,  (ID 


where 


•RVB*  "  (R/V)K  “  co*  •  «<•»  *>  (12) 

g(e*  4)  •  £(sin  e  sin  #)2  +  cos2  ej  (13) 

As  indicated  by  Alpers  and  Rufenach  (1979),  the  effect  of  nonlinearity  can  be 
to  enhance  periodic  image  intensity  (i.e.,  accentuate  wavelengths)  as  shown  in 
figure  8,  but  extreme  care  must  be  taken  since  this  accentuation  will  be 
displaced  in  azimuth  and  phase. 

For  azimuth- traveling  waves  (4*0)  and  the  previously  selected  examples, 
the  maximum  peak- to- trough  modulation  is  estimated  (Ao  -  <o>  1  +  C)  as  11  dB 
for  SEASAT  (X  -  100  m,  A  -  0.5  m)  and  7.1  dB  for  JOWIP  (X  -  35  mT  A  -  0. 1  m) . 
The  corresponding  values  of  the  linearity  parameter  are  C  ■  2.7  and  C  -  1.1, 
respectively.  Figure  9  summarizes  the  relative  variation  of  three  image- 
modulation  transfer  functions  with  long  wave  propagation  direction.  Table  4 
assembles  the  maximum  crest- to- trough  RCS  values  calculated  above  for  easy 
comparison.  Note  that  these  maximum  values  were  calculated  using  absolute 
values  of  the  complex  FFTs;  they  cannot  be  easily  combined. 
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to  resolve  at  25  meters  although  line  enhancement  might  adequately  define 
elongated  wave  crests . 


Figure  7.  Modulation  transfer  function 
R  ^  for  the  tilt  mechanism  as  a  func¬ 
tion  of  incidence  angle  for  HH  and  W 
polarization.  The  upper  curves  <  —  R ||  /IK  || 
apply  to  range- traveling  waves  and  the 
lower  curves  are  for  azimuth-traveling 
waves.  These  calculations  apply  to  an  L- 
band  radar  as  used  for  SEASAT  SAR.  The 
SEASAT  SAR  angle  of  incidence  8  ■  23* 

( adapted  from  Alpers  et  al . ,  1981). 


INCIDENCE  ANGLE  (0,1 


For  the  JOWIP  radar,  the  estimated  tilt  modulation  for  an  A  -  0.1  m,  K  - 
2v/35  -  0.18  range- traveling  wave  would  be  to  •  <oQ>  [1  +_  1.8  (0.1)],  which 
corresponds  to  a  maximum  tilt  modulation  of  1.4  dB  crest-to- trough. 

For  azimuth- traveling  ships,  the  estimated  modulations  are  too  small  to 
be  important;  0.2  dB  and  0.1  dB  for  the  SEASAT  and  JOWIP  cases,  respectively, 
which  is  below  the  expected  coherent  speckle  level  (~  0.5  dB).  For  a  ship 
heading  at  an  angle  of  55  degrees  relative  to  the  SAR  track,  the  range 
traveling  side  of  the  wake  envelope  (a  -  35*)  should  display  a  much  sharper 
tilt  modulation  than  the  opposite  wake  ora  (Kelvin  envelope  propagation 
direction  -  ship  heading  +_  35*). 

HYDRODYNAMIC  MODULATION 

In  the  first  order  theory  for  decimeter  waves  (Alpers  et  al.,  1981), 
modulation  results  from  purely  hydrodynamic  interactions  between  the  short 
Bragg  scattering  waves  and  the  longer  ocean  wavelengths.  Nonuniform  wave- 
induced  airflow  is  neglected.  When  the  spectrin  of  short  surface  waves  is 
assumed  to  be  a  Phillips  spectrum  ( E  ~  |  K| ~4 ) ,  and  the  Bragg  resonant  wave¬ 
length  (fc)  is  taken  to  be  in  the  range  direction,  the  MTF  can  be  written  in 
terms  of  the  angle  (•)  between  the  long  wave  (K)  and  the  SAR  ground  track 
(Alpers  et  al.,  1981), 


14 


o  -  0  ♦  -  o0  [l  ♦  Jl R(K)  Z(K) 


i(KX  -  Mt) 


♦  CC]  «J  , 


where 


R(K)  -  R(K)  ♦  R(K)_  .  ♦  R(K)„  ,  .  k,  .  <«) 

tilt  Hydro  Vel  Bunching 

The  complex  Modulation  transfer  function  (MTF) ,  R(K),  includes  the  combined 
effects  of  the  three  mechanisms  Mentioned  shove.  Hare  is  the  cross-section 
associated  with  local  surface  roughness;  Z(K)  is  the  Fourier  transfore  of  the 
amplitude  (A)  of  the  ocean  surface;  K  and  a  are  wove  vectors  and  radian  fre¬ 
quencies  of  the  large  scale  wave  field;  and  CC  is  the  coaiplex  conjugate. 

If  only  one  long  wave  (K)  is  presented  and  the  linearity  conditions  to  be 
discussed  below  are  net,  then  this  expression  can  be  simplified  to  read 


1  +  R_ 


Ar  cos  (Kx  -  ut  +  4)j 


6  -  tan 


■1  [In  Rl 

Lrs  rJ 


irized  below 


RK  "  Rtilt  +  RHydro  +  RVel  Bunching, 

facilitating  meaningful  comparison  of  the  individual  MTF's 
for  specific  wave  directions  and  amplitudes. 

TILT  MODULATION 


Assuming  a  uniformly  wind-roughened  surface,  the  effect  of  slope  varia¬ 
tions  at  different  locations  along  the  long  wave  gives  rise  to  the  theoretical 
MTF  estimates  in  figure  7  (Vesecky  and  Stewart  (1982)).  The  two  parts  of 
Rtllt>  R  ||  and  Rj^,  refer  to  tilt  components  in  and  perpendicular  to  the  plane 
of  incidence.  For  SEASAT/JOWIP  incidence  of  23  degrees,  we  have  for 
range- traveling  waves  (•  -  0*),  R  *  R||  *  1°  K;  and  for  azimuth-traveling 
waves  R  *  Rx  ^  0.6  K.  The  maximum  modulation  occurs  for  waves  propagating 
in  the  range  direction. 

For  SEASAT,  the  estimated  maximum  tilt  cross-section  modulation  for  an  A 

-  0.5  meter,  K  -  2w/100  »  0.06  range- traveling  wave  would  be  Ac  -  <oQ>  [1  ± 

R^ilt  A]  -  <oQ>  [1  ±  0.6  (0.5)],  which  corresponds  to  2.4  dB  crest- to- trough 

ratio.  This  level  of  modulation,  possibly  indicative  of  the  transverse  waves 
behind  a  large,  high-speed  (25  knot),  range-traveling  ship  would  be  difficult 


SAR  IMAGE  MODULATION 


For  the  low  to  moderate  surface  wave  conditions  in  the  JOWIP  test  areas, 
the  linear  theory  summarized  by  Alpers  et  al.,  (1981)  is  the  only  comprehen¬ 
sive  theory  available  for  describing  SAR  ocean  wave  image  modulation 
mechanisms.  This  theory  is  based  on  a  a  two- scale  model  of  the  ocean  surface 
(Wright,  1968).  Using  this  approach,  longer  gravity  waves  modulate  the  micro¬ 
wave  backscatter  due  solely  to  local  wind-generated  surface  roughness  which  is 
assumed  to  be  nearly  isotropic  and  can  be  described  in  terms  of  a  saturated 
wave  number  spectrum  (Phillips,  1977).  This  roughness  is  attributed  specific¬ 
ally  to  those  short  gravity  waves  traveling  radially  with  respect  to  the  radar 
and  obeying  the  Bragg  resonance  condition  K  -  2K  sin  0;  where  K  is  the  ocean 

wave  number,  K  is  the  radar  wave  number,  and  8  is  the  local  incidence  angle. 

R 

The  arousal  time  for  SEASAT  Bragg  resonant  waves  (  ~30  cm  @  23*)  is 
estimated  to  be  quite  short:  approximately  25  seconds  for  an  energy  increased 
by  a  factor  of  3  for  a  5  m/s  wind.  If  the  wind  speed  is  referenced  to  10 
meters  above  the  surface  and  dissipative  forces  are  ignored  (Vesccky  and 
Stewart,  1982),  this  time  varies  inversely  with  the  square  of  wind  speed 
between  2  and  20  m/s . 

Dissipation  of  these  30-cm  waves  is  poorly  understood  and  depends  on 
wave-wave  interactions,  viscous  damping,  and  breaking  to  sustain  wave  number 
saturation  of  the  power  spectrum  under  conditions  of  increased  wind  energy 
input.  In  the  absence  of  wind  input  and  wave-wave  interaction,  the  clean 
water  viscous  energy  density  decay  time  (e- folding)  for  SEASAT  resonant  (30 
cm)  waves  has  been  roughly  estimated  as  450  seconds.  In  the  real  ocean  this 
will  vary  widely  and  could  be  an  order  of  magnitude  shorter  with  the  intro¬ 
duction  of  surface  films  (Phillips,  1977)  or  strongly  interactive  waves. 

Although  alternative  formulations  for  time  variant  SAR  scenes  are  possi¬ 
ble  (Hager,  1980;  Shuchman  and  Zalenka,  1978;  and  Shuchman,  1981),  the 
specific  long  ocean  wave  image  modulation  mechanisms  supported  by  general 
consensus  in  the  radio-oceanographic  community  (Vesecky  and  Stewart,  1982) 
are: 

e  tilt  of  the  scattering  surface  (Valenzuela,  1978); 

e  hydrodynami ca 1 ly  induced  crest-to-trough  variations  in  the  density  of 

present  scatterers  (Keller  and  Wright,  1975);  and 

e  apparent  crest  and  trough  coordinate  distortions  (velocity  bunching) 
caused  by  long  wave  orbital  velocities  (Alpers  et  al.,  1981). 

The  latter  mechanism  can  result  in  contrast  enhancement  for  long  wave  trains 
but  could  also  smear  out  the  relatively  lonesome  wave  components  of  a  ship 
wake  (Alpers,  1983). 

The  normalized  radar  cross-section  (o)  modulation  by  long  ocean  waves  is 
statistically  described  as 
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Imagt  Friction  (K) 


Figure  5.  Ground  incidence  angle  (8^) 
and  ground  range  (D)  aa  a 
function  of  fractional  dis¬ 
tance  (K)  froa  near  edge  on 
slant  iaages. 
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Figure  6.  Slant-to-ground  angle 
conversion  for  JOWXP 
iaages  as  a  function  of 
surface  incidence  angle 


The  imaging  geometry  for  the  CCRS/ERIM  airborne  system  is  shown  in  figure 
4.  Since  the  primary  radar  coordinates  are  range  and  azimuth  angle,  all 
processing  is  referred  to  the  slant  or  range  azimuth  plane  as  shown  in  figure 
4  where  the  range  delay  and  5. 4-tan  sweep  translate  into  a  distorted  coverage 
over  approximately  10  km  on  the  ground.  At  the  -degree  nominal  incidence 
used  for  comparisons  with  SEASAT,  the  range  distances  are  considerably  fore¬ 
shortened  (K  »  0. 1).  These  foreshortened  distances  yield  small  slant  plane 
angles  (3  degrees  vs  8  degrees)  which  must  be  corrected  to  the  ground  plane 
for  comparison  with  the  familiar  SEASAT  surface  plane  (4-look)  nominal  25 
meter  x  25  meter  resolution  images.  Since  all  but  one  of  the  JOWIP  images 
presented  here  will  be  slant-plane ,  figures  5  and  6  are  included  to  facilitate 
comparisons . 

Figure  5  relates  the  fractional  distances  (K)  (from  the  near  edge  of  the 
image  to  the  target  area)  to  nominal  local  incidence  angle  (8^  and  ground 
distance  from  nadir  (D) .  Figure  6  can  be  used  for  conversion  of  slant  to 
ground  angles . 


Platform  JOWIP  flight/SAR  parameters 


0.0  Image  fraction  IK)  1  .00 

0.7S  km  Surface  distance  (D)  10.3  km 

6.1°  Incidence  angle  ($ | )  55.7° 


Figure  4.  JOWIP  surface  incidence  angles  and  distances. 
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Both  the  phase  velocity  and  N(z)  increase  as  the  density  gradient  dp/dz  in¬ 
creases.  In  the  case  of  the  QOAPAW,  intended  ship  speeds  were  on  the  order  of 
14  knots  (7  b/s)  which  is  such  greater  than  possible  phase  velocities  asso¬ 
ciated  with  wavelengths  near  the  ship  length  of  60  n.  Here  ship  speed  is 
above  the  phase  velocity  and  a  narrow  vee  pattern  would  be  produced  behind  the 
ship.  The  smaller  the  phase  velocity  (lower  density  gradient),  the  smaller 
the  angle  produced  in  the  vee  pattern. 

Mien  the  turbulent  wake  is  formed,  the  existing  density  gradient  is 
temporarily  mixed.  The  density  of  wake  water  near  the  top  of  the  wake  is 
greater  than  that  outside  the  wake  and  the  reverse  is  correspondingly  true  at 
the  bottom  of  the  wake.  As  this  water  returns  to  equilibrium,  internal  waves 
are  radiated.  A  vee  pattern  is  again  produced  by  waves  traveling  away  from 
the  ship's  speed.  The  amplitude  of  these  waves  depends  on  the  magnitude  of 
the  density  gradient  across  the  wake  and  the  final  height  of  the  msks  before 
collapse.  For  expected  Da bob  Bay  density  gradients  the  vee  angle  pattern  for 
QOAPAW  was  predicted  to  be  1  degree  or  less  with  vanishingly  small  amplitude. 

SAR  PARAMETERS  AND  IMAGE  FORMAT 

Since  the  primary  objective  of  the  ship  wake  investigation  was  to  explain 
SEASAT  images  it  was  necessary  to  adhere  as  closely  as  possible  to  the  SEASAT 
imaging  geometry.  This  involved  many  tradeoffs  including  consideration  of 
aircraft  capabilities,  surface  radar  incidence,  small  signal  suppression,  and 
emphasis  on  the  horizontally  polarized  L-band.  The  primary  flight  and  radar 
parameters  for  the  SEASAT  and  OCRS  SAR  ship  wake  images  were  selected  by  NOSC 
and  ERIM  and  are  shown  in  table  3  for  comparison  purposes. 


Table  3.  Comparative 

SAR  image 

parameters 

• 

PARAMETER 

SEASAT 

CCRS  SAR 

SAR  wavelength  (L-band)  (m)  XQ 

0.235 

0.235 

Ground  velocity  (km/s)  V 

7.2 

0.13 

Altitude  (km)  H 

782 

7.01 

Surface  incidence  (deg)  4 

(ship  wakes  usually  generated  at 
nominal  SEASAT  23*  incidence) 

19  - 

25 

6-55 

Mean  slant  range  R  “  H/cos  23  (km) 

850 

7.6 

Azimuth  beamwidth  (deg)  0 

A 

1 

6 

Bandwidth  (MHz)  Br 

19 

60 

Ground  beamwidth  (Ion)  L  «  R0a 

89.0 

(14  used) 

0.796 

Azimuth  resolution  (m)  (1-look) 
pa  -  XQ  R/(2TV) 

6.9 

5.0 

Range  resolution  (m)  pr  *  C/(2Br) 

7.9 

3.0 

Coherent  integration  time  (s) 

T  -  XQ  R/(2paV) 

2.0 

1.4 

B.  JOWIP  SAR  Plight  9,  Pass  1,  28  July  1983,  13:17:08  POT 


Figure  12.  Range- traveling  ship  wake  comparison. 
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estimated  by  densitometer  scans  across  the  SAR  image  to  be  10  dB  above  the 
local  background  (which  is  estimated  at  3  dB  above  system  noise).  Also  note 
that  the  wake  shows  periodic  modulation  which  closely  approximates  (40  m)  the 
transverse  wavelength  associated  with  ship  speed  (15  knots),  dose  inspection 
reveals  that  the  port  and  starboard  modulations  are  out  of  phase  revealing  a 
measurable  relative  azimuth  displacement. 

Figure  11  shows  similar  but  smaller,  narrow-vee,  full  angles  of  5.0  and 
3.6  degrees  for  the  SEASAT  (23°  incidence)  and  JGMIP  (50°  incidence)  SARs, 
respectively.  These  ship  headinqs  make  an  angle  of  nominally  30  degrees  with 
the  SAR  ground  track.  Note  that  the  50-degree  incidence  for  the  bottom  image 
makes  specular  scattering  extremely  unlikely  as  well  as  indicating  less  di¬ 
rectionality  for  the  ship-induced  roughness.  While  the  lower  wake  image  shows 
a  measured  modulation  length  of  35  meters  (tfaich  could  correspond  to  the 
transverse  waves  associated  with  a  14-knot  ship),  the  upper  SEASAT  wake  image 
modulation  length  is  a  long  140  meters.  This  case  matches  more  closely  the 
local  swell  wavelength  measured  in  the  wind-roughened  area  in  the  upper  left 
and  propagating  at  a  si  35  degrees  relative  to  the  SAR  ground  track. 

R>r  the  range-traveling  ship  comparisons  shown  in  figure  1 2,  the  wake 
arms  are  completely  absent  and  the  only  remaining  wake  characteristic  is  the 
turbulent  wake  aft  of  the  ship.  Ebr  the  top  SEASAT  image  this  appears  as  the 
usual  dark  (lack  of  return)  against  a  wind- roughened  background.  In  the  lower 
J0WIP  image,  the  turbulent  wake  appears  as  a  bright  uniform  return  across  an 
approximately  constant-width  turbulent  wake  viewed  against  the  dark  background 
for  this  calm  wind  (  ~0.5  m/s)  case.  The  sharp  returns  defining  the  usual 
wake  arms  are  absent  with  reduced  image  density  ~5  dB  above  background. 
Close  inspection  of  this  wake  failed  to  reveal  a  definite  wake  modulation 
length  in  the  bright  ship-roughened  area.  Table  8  summarizes  the  JOWIP  ship 
wake  measurements  by  date,  flight,  and  pass  number.  The  ship  headings 
( |  Hq|  <1 80° )  are  measured  from  the  SAR  track  to  the  ship  track  and  labeled 
negative  (-)  if  closing  and  positive  (+)  if  opening  range.  Angles  are 
positive  counterclockwise  for  the  right-looking  antenna  and  positive  clockwise 
for  the  left-looking  antenna.  The  not  resolvable  (NR)  wake  characteristics 
will  all  be  roughly  explained  in  terms  of  ambient  wind  and  ship  heading. 

All  the  SAR  passes  except  1-4,  7-1,  7-3,  and  9-1  were  intended  for 
azimuth-traveling  ship  wakes.  Run  1-4  resulted  from  an  unknown  volunteer  ship 
and  the  latter  passes  were  intended  to  image  range- traveling  ship  wakes. 

In  only  one  case  in  a  localized  transitional  wind-roughened  area  during 
relatively  calm  conditions  (averaged  2.5  m/s)  did  the  L-band  SAR  produce  a 
recognizable  Kelvin  wake  pattern.  This  pattern  is  shown  in  figure  13  where  it 
is  compared  with  the  only  known  SEASAT  image  displaying  the  transverse  wave 
pattern  within  the  Kelvin  wake  envelope  angle. 

The  SEASAT  image,  processed  by  the  Royal  Air  Establishment,  to  21  by  12 
meter  range  and  azimuth  resolution,  displays  a  94-meter,  approximately  range¬ 
traveling  (♦  *  1 04# ) ,  transverse  wavelength  which  compares  favorably  with  the 
measured  azimuth  displacement  (&x  »  570  m)  for  a  23-knot  ship*  v  -  Ax  V/(  R 
sin  4  sin  9).  Even  with  the  vast  ocean  surface  coverage  attained  by  the 
SEASAT  SAR,  high  speed  targets  such  as  this ,  moving  in  the  range  direction  to 
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SEAS AT  SAR  Rev  834,  24  Aug  1978,  07:27:12  GMT 
(Courts ay  of  RAE  Fambourough,  England) 

SAR  TRACK 


B.  JOW2P  SAR  Plight  8,  Pass  4,  28  July  1983,  09:47:07  PDT 
Figure  13.  SEASAT  and  JOWZP  transverse  wavs  pattern  conparison. 
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minimize  wave  motion  near  effects  (Alpers  et  al.,  1981),  are  a  rare  com¬ 
modity.  The  average  wind  speed  for  this  image  is  unknown. 

The  CCRS/ERIM  image  on  the  lower  half  of  figure  13  also  displays  the 
transverse  wave  pattern  with  wake  half  angles  (  —19°)  nearer  the  theoretical 
Kelvin  envelope.  However  this  image  is  for  an  azimuth-traveling  ship  passing 
through  slightly  wind- roughened  water  with  a  surface  radar  incidence  angle  of 
approximately  41  degrees,  a  relative  ship  heading  of  19  degrees,  and  a 
measured  transverse  wavelength  of  39  meters.  Here  the  entire  transverse  wave 
pattern,  slightly  limited  in  spots  by  insufficient  surface  roughness,  is 
visible . 


DATA  ANALYSIS  AND  DISCUSSION 


This  section  covers  special  analysis  efforts  of  8AR  iaage  data,  oblique 
photography,  vertical  photography,  stereo  mensuration,  inner  sake  slope 
observations,  and  an  analysis  smsary. 

SAR  IMAGE  DATA 

Reordering  the  SAR  wake  measurements  presented  in  table  8  by  the  local 
wind  speeds  shown  in  table  6  gives  us  table  9  which  clearly  shows  the  effects 
of  local  wind  speed  on  JOWIP  SAR  wake  component  visibility. 

Here  wind  speeds  have  been  grouped  into  three  categories  according  to 
wake  component  visibility.  For  the  first  essentially  flat  cals  category,  the 
dark  (lack  of  return)  turbulent  wake  is  not  seen.  On  run  9-1  the  turbulent 
wake  is  clearly  due  to  the  increased  scatter  generated  by  the  range-traveling 
ship  as  are  the  arms  apparent  for  the  asiauth-traveling  ship  on  run  9-2. 

fbr  the  transitional  wind  category,  the  situation  is  sosewhat  sore  com¬ 
plex  with  turbulent  wake  measurement  complicated  by  the  ship  passing  in  and 
out  of  small  local  wind-roughened  patches.  The  important  observation  here  is 
that  under  these  wind  condition*  the  ship-generated  surface  roughness  persists 
in  the  locally  calm  areas  along  the  track  and  gives  rise  to  measurable  narrow 
wake  angles,  These  narrow  angles  are  further  seen  to  exist  for  relative  ship 
heading  out  to  at  least  35  degrees  and  for  radar  surface  incidence  angles 
between  23  and  50  degrees.  In  only  one  case  (run  8-4)  at  the  upper  end  of 
this  wind  category,  was  a  wind- gene  rated  surface  roughness  of  2.5  m/s  barely 
adequate  to  reveal  the  full  width  of  the  transverse  pattern  between  the  Kelvin 
wake  envelope  arms.  The  wind  speed  of  3.5  m/s  is  the  estimated  speed  required 
to  maintain  a  saturated  a pec trim  of  first  order  Bragg  scattering  (30  cm) 
gravity  wave  roughness  (Vesecky  and  Stewart,  1982). 

The  transverse  ship  wake  waves  in  the  wind-roughened  category  were 
probably  not  of  sufficient  amplitude  (£  0.5  m)  to  be  visible  by  the  usual 
iaage  modulation  mechanisms  (Vesecky  et  al.,  1982).  The  last  trace  of  the 
wake,  the  turbulent  patch,  disappears  at  the  higher  wind  speeds  (>10  m/s) 
making  even  the  SAR  determination  of  ship  heading  impossible. 

SAR  WAKE  IMAGE  MODULATION  LEVELS 

In  order  to  quantify  the  observed  iaage  modulation  levels,  a  series  of 
cross  wake  (radar  range  direction)  scans  were  run  on  an  ERIM  optical  pro¬ 
cessor.  Ptor  this  evaluation  a  calm  surface  portion  of  the  wake  on  SAR  run  8-2 
between  800  and  1200  meters  aft  of  QUAPAW  was  selected.  The  scans  were  accom¬ 
plished  in  the  slant  range  plane  using  an  8-aeter-diameter  aperture;  each  line 
scan  was  plotted  as  linear  intensity  versus  slant  range  to  explore  the  effects 
of  observed  transverse  wave  (1  -  40  m)  modulation  over  a  10-wavelength  inter¬ 
val  (51  scan  lines).  Figure  14  shows  an  example  of  five  of  the  51  consecutive 
scans  across  the  wake  with  an  approximately  8-meter  along-track  spacing 
between  scans.  The  side-to-side  asymmetry  is  apparent  as  are  multiple  smaller 
peaks  in  addition  to  the  predominant  ones  responsible  for  narrow  wake  arms. 
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Figure  14.  Linear  intensity  versus  radar  range  across  QUAPAW  wake 
(Run  8-2);  ERIM  optical  processor,  26  March  1984. 


Statistics  for  the  entire  51 -scan  data  set  are  shown  schematically  in 
figure  15.  In  treating  these  data  the  wake  centerline  was  determined  for  each 
scan  line  by  bisecting  the  distance  between  the  major  peaks.  The  distances  of 
each  of  the  peaks  from  this  centerline  were  then  measured  and  they  consis¬ 
tently  fit  into  reasonably  well  defined  groupings  centered  about  the  wake 
coordinate  angles  shown  with  the  frequencies  indicated.  The  amplitudes  ware 
then  measured  for  the  background  and  each  peak  along  with  these  linear  inten¬ 
sity  values  converted  to  the  mean  decible  values  shown  in  the  figure.  Note 
that  a  peak  (AV  6.4  dB)  occurred  between  the  major  peaks  in  57  percent  of  the 
cases  and  that  in  the  remaining  43  percent  of  the  scans  the  centerline  image 
intensity  was  3.7  dB  above  this  locally  calm  background  level.  Me  suggest 
that  a  wind-induced  7-dB  increase  in  the  background  level  would  eliminate  all 
but  the  major  peaks  and  result  in  the  usual  dark  (lack  of  return)  turbulent 
wake  along  the  ship's  track.  Of  possible  significance  are  the  very  low  fre¬ 
quency  and  amplitude  associated  with  the  Kelvin  wake  envelope  and  the 
generally  higher  amplitudes  on  the  far  range  side  of  the  wake. 


Port  Arm 
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Figure  15.  SAR  wake  image  scan  statistics  for  calm  water  (SAR  pass  8-2). 


Along- track  spacing  of  peaks  was  estimated  from  along-ship  track  wake-arm 
fluctuations  in  image  intensity  to  be  (43  port  +  46  starboard) /2  ■  44.5 
meters.  Attempts  to  quantify  the  side-to-side  phase  difference  have  been  less 
successful  but  something  close  to  half  the  transverse  wavelength  appears 
probable  from  inspection  of  images.  If  the  ship- induced  roughness  is  asso¬ 
ciated  with  short  gravity  waves  propagating  in  roughly  opposite  crosstract 
directions,  then  a  phenomenon  analogous  to  doppler-splitting  (Alpers  et  al., 
1981)  might  cause  a  nominal  15-meter  azimuth  offset  of  the  amplitude  peaks. 

It  must  be  pointed  out  that  the  along-track  wake  arm  modulations  varied 
between  7  and  16  dB  on  the  far  range  side  of  the  wake  and  between  6  and  14  dB 
on  the  near  range  side.  These  peak- to- trough  differences  of  9  and  8  dB  are 
far  in  excess  of  the  theoretical  predictions  for  tilt  or  hydrodynamic  image 
modulation  but  very  close  to  the  velocity-bunching  modulation  level  (7  dB) 
predicted  in  table  4.  We  therefore  conclude  that  an  effect  like  "velocity 
bunching,"  as  described  by  Alpers  (1979,  1981,  1983),  caused  by  transverse 
wave  orbital  notion,  is  the  SAR  image  modulation  mechanism  responsible  for 
narrow  wake  images. 

OBLIQUE  PHOTOGRAPHIC  OBSERVATIONS  OP  SURFACE  ROUGHNESS 

Figures  16  and  17  illustrate  the  effects  of  low  wind  speed  surface  rough¬ 
ness  on  a  SAR  image.  The  SAR  image  in  figure  16A  covers  a  7. 4- kilometer. 
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Turbulent  Wake 


Figure  17B.  Photo  of  surface  roughness  in  the  vicinity  of 

CPAV  ENDEAVOR.  DEEP  PHOTO:  728081848 


41 


left-looking  slant  range  segment  of  run  8-2.  The  image  times  were  approxi¬ 
mately  0815:30  to  0816:30  PDT  and  the  average  wind  speed  recorded  aboard 
ENDEAVOUR  (shown  in  figures  16A,  16B,  and  17B)  over  the  interval  0812-0817  was 
1.5  +  0.3  m/s,  temporal  average.  The  SAR  image  suggests  a  considerable 
spatial  variation  in  surface  roughness  and  the  intent  of  this  example  was  to 
correlate  qualitative  estimates  of  surface  roughness  from  photographic 
coverage  with  the  SAR  image. 

Six  aerial  photos  were  shot  from  1000  feet,  over  an  interval  of  5.5 
minutes  spaced  around  the  1 -minute  necessary  for  the  SAR  image.  Figure  16B 
shows  a  composite  of  subjective  surface  estimates  on  a  partial  chart  of  Dabob 
Bay.  The  wind-roughened  areas  were  estimated  from  six  oblique  aerial  photos; 
the  two  most  helpful  are  shown  as  parts  A  and  B  of  figure  17.  Figure  17  shows 
the  scar  caused  by  passage  of  the  QUAPAW  through  a  locally  wind-roughened  area 
looking  southeast  1 -minute  before  SAR  time.  Notice  how  the  turbulent  wake  has 
calmed  the  wind-roughened  area  in  the  upper  left  corner;  notice  also  the  sharp 
edge  between  the  smooth  and  roughed-up  areas.  This  transition  between 
roughened  and  calm  areas  is  remarkably  abrupt  and  these  abrupt  lines  strongly 
persist  from  photo  to  photo  over  many  minutes  of  elapsed  time,  resulting  in  a 
quasi-standing  rather  than  a  propagating  surface  wave  disturbance.  The  sharp 
edges  of  these  patterns  are  easily  seen  on  the  SAR  image  as  it  slants  off  to 
the  upper  left  corner.  Figure  1 7B  was  imaged  over  2  minutes  after  the  SAR 
pass,  and  shows  water  surface  detail  near  the  Canadian  ship  ENDEAVOUR  as  it 
moves  slowly  across  the  QUAPAW's  wake  at  a  heading  of  090°T.  Notice  how  the 
composite  diagram  of  roughened  surface  areas  from  the  oblique  photos  closely 
matches  the  SAR  image  pattern  of  the  same  scene.  Notice  also  that  since  the 
turbulent  wake  strip  produces  a  lack  of  return  on  the  SAR  image,  it  is  not 
visible  in  the  calm  area  between  two  narrow-angle  arms  of  the  wake,  but  is 
only  evident  when  there  is  a  sufficiently  roughened  surface  to  show  this  lack 
of  radar  return.  These  roughened  background  areas  are  thought  to  be  populated 
by  capillary  and  short  gravity  waves  with  centimeters  amplitudes. 

VERTICAL  PHOTOGRAPHY  FOR  SHORT  GRAVITY  WAVELENGTH  AND  HEIGHT  DETERMINATION 

Vertical  photographs  taken  over  the  Dabob  Bay  SAR  events  fall  into  two 
altitude  regimes.  We  took  a  sequence  of  color  runs  covering  the  entire  wake 
system  in  a  series  of  single  frames  (135  m  x  135  m)  taken  from  a  750-foot 
altitude,  only  during  DREP  run  9-1.  The  bulk  of  the  photography,  however,  was 
high  resolution  black  and  white  stereo  pairs  covering  small  sectors  (60  m  x  60 
m)  of  the  wake  system  taken  from  a  200-foot  altitude.  All  frames  are  nomi¬ 
nally  2x2  inch  format  (70  mm  film).  Using  both  these  phoco  sources  as 
single  frames  for  wavelength  extraction  will  be  discussed  here  and  a  brief 
description  of  progress  on  stereo  analysis  of  instantaneous  ocean  surface 
topography  will  follow  in  the  next  section. 

One  of  the  color  photo  runs,  a  short  sequence  of  frames  on  both  the  port 
and  starboard  sides  of  the  wake  pattern  approximately  1000  meters  aft  of  the 
ship  was  used  for  preliminary  estimate  of  wake  wavelengths  and  propagation 
directions  with  respect  to  the  ship's  track.  These  parameters  are  plotted  in 
figure  18  where  they  are  compared  with  the  calculated  Kelvin  wake  parameters 
defined  by  equations  (2)  and  (3).  Wake  wavelengths  about  10  meters  on  this 
calm  wind  interval  compare  favorably  with  the  calculated  Kelvin  wake  values, 
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Figure  18.  Preliminary  photographic  estimates  of  ship  wake  component 
wavelength  vs  relative  propagation  angle  (run  9-1). 


but  below  10  meters  all  measurable  wavelengths  were  equal  to  or  greater  than 
the  calculated  Kelvin  theory  predictions. 

Also  included  in  figure  18  are  several  ambient  wave  components  which 
compare  favorably  with  the  average  wind  direction,  also  indicated  in  the 
figure.  This  may  be  fortuitous  at  these  low  wind  speeds  (<  1  m/s)  where  wind 
direction  commonly  wanders . 

The  observed  surface  wavelength  deviations  from  the  Kelvin  theory  could 
result  from  (1)  background  wave  clutter;  (2)  inappropriate  directional  sun-sky 
illumination  at  various  phases  (slopes  along  the  longer  transverse  wave¬ 
lengths);  (3)  reduced  slopes,  and  therefore  reflected  light  from  these  shorter 
wavelengths;  or  (4)  failure  of  simple  Kelvin  model  predictions  at  these  short 
wavelengths . 

It  must  be  emphasized  that,  while  photographically  observed  wavelengths 
and  propagation  directions  are  valid,  the  same  wavelength  propagating  in 
another  direction  night  not  be  as  obvious.  Further,  we  cannot  state  that 
because  shorter  wavelengths  are  not  apparent  on  the  photos,  they  are  not 
present  in  the  scene. 
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Since  items  (1)  and  (2)  above  ware  under  reasonable  control  by  site 
selection,  ship  track  orientation,  and  luck,  items  (3)  and  (4)  (which  must 
also  be  considered  for  all  the  photographic  records)  were  felt  to  be  the 
critical  factors.  Approximate  sun  azimuth  and  elevation  for  each  QUAPAW  wake 
run  are  given  in  table  10.  For  vertical  photography,  the  table  also  gives  the 
critical  slope  angles  required  for  direct  sun  glint  or  glitter,  equal  to  half 
the  sun's  zenith  angle.  Problems  associated  with  using  the  width  of  the 
glitter  patterns,  which  vary  with  surface  slope  (Hulbert,  1934)  and  with  the 
statistical  extraction  of  two-dimensional  slope  spectra  from  within  glint 
pattern  photographs  (Cox  and  Munk,  1954),  were  to  be  avoided  here  by  using  a 
method  analogous  to  the  optical  transform  method  of  Stillwell  (1969).  By 
working  outside  the  glitter  pattern  and  relying  on  a  continuous  skylight 
luminance,  we  hoped  that  surface  measurements  could  be  made  with  a  minimum  of 
directional,  slope,  and  wavelength  bias. 


Table  10.  Sun  azimuth  and  elevation  for  photographic 
ship  wake  runs  (47-41N;  122-1 5W) . 


Run 

Time  (PDT) 

Sun  Azimuth 
(deg) 

Sun  Elevation 
(deg) 

Critical  Angle 
(deg) 

7-1 

16:18:54 

247 

46 

22 

8-1 

07:30:44 

79 

16 

37 

8-2 

08: 16:09 

88 

24 

33 

8-3 

09:01:18 

96 

31 

30 

8-4 

09:47:04 

105 

39 

26 

9-1 

13:16:13 

178 

61 

15 

9-2 

14:00:48 

200 

60 

15 

9-3 

14:46:22 

218 

56 

17 

Using  the  sequence  of  vertical  photos  from  via  low  altitude  run  (8-3) 
positioned  in  wake  coordinates  in  figure  19,  the  subdivided  (10  m  x  10  m) 
regions  between  grid  marks  were  digitized  using  a  50-micron  aperture  scanning 
densitometer  on  the  70-mm  original  negatives.  From  this  process  digital  two- 
dimensional  fast  Fourier  transforms  ( FFT)  of  the  192  x  192  square  array  were 
obtained.  The  associated  FFT  numbers  are  also  shown  in  the  figure  along  with 
a  typical  frame  ( •  1 3 ) .  The  entire  run,  frames  1-23  covered  the  path  shown 
from  900  meters  aft  to  over  300  meters  ahead  of  the  snip.  Because  the  wake 
center-line  was  visible  and  the  ship  travel  speed  and  aircraft  flight  param¬ 
eters  known,  this  photo  track  was  established  at  a  10-degree  crossing  angle. 
Adjusting  this  flight  line  from  ground  coordinates  to  the  stationary  ship  wake 
coordinate  system  changes  this  crossing  angle  to  the  1 3  degrees  shewn  on  the 
drawing.  The  wake  angles  and  distance  aft  of  the  ship  corresponding  to  each 
FFT  are  included  in  the  data  tables  to  follow. 
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Table  15. 

Estimated  surface  wavelengths  for  JONIP  inner 

wake  arms 

First  order  Bragg 

Nearest  cross¬ 

SAR  Run 

wavelength  (cm) 

track  (cm) 

Other  arm 

8-1* 

24 

50+20 

31+20 

8-2 

30 

63^30 

63+30 

8-3 

30 

4100 

41+30 

8-4 

18 

95+30 

44+30 

9-2 

31 

44+30 

44+30 

9-3 

26 

44+30 

41+30 

*V  «  6.5  m/s 


^SAR  TRACK 


Figure  27.  BRIM  digitally  processed  JCWIP  ship  wake  linage  (run  8-2). 


should  he  tested  against  other  extreme  parametric  variations.  Figures  28A  and 
B  show  a  slant  range  comparison  of  a  narrow  wake  image  made  by  the  same  BRIM 
SAR  in  1978  with  both  X-  and  L-bands  functioning  properly.  While  the  L-band 
shows  a  much  stronger  wake  with  a  4.7-degree  (ground  range)  wake  angle,  at 
this  steep  incidence  (18°)  the  X-band  also  shows  a  shorter  narrow  wake  at  3.7 
degrees  (ground  range  half  angle),  which  is  modulated  by  a  54-meter  transverse 
wavelength  (17.8-knot  ship).  For  x-band  at  this  incidence,  the  Bragg  length 
is  approximately  5  cm  and  the  wake  extends  for  2.5  km. 
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Table  14.  Comparison  of  wake  angles  (deg)  with  Bragg  resonant 
wavelengths  (cm)  for  various  incidence  angles  (deg) 


(Azimuth  Traveling  Ships  4  -  o) 


Radar  Incident  Angle  (8^) 

20  23 

25 

30 

40 

50 

90 

1st  order  Bragg  wavelength  (1Q) 

34.4  30.1 

27.8 

13.5 

18.3 

15.3 

11.8 

Kelvin  wake  angle  (8R) 

2.6  2.5 

2.4 

2.2 

1.9 

1.7 

1.5 

Cross- track  roughness  angle  (6  ) 

2.7  2.5 

2.4 

2.2 

1.9 

1.8 

1.6 

2nd  order  Bragg  wavelength  (X  ) 

2o 

68.8  60.2 

55.4 

47.0 

36.6 

30.6 

23.6 

Kelvin  wake  angle  (6g) 

3.7  3.5 

3.4 

3.1 

2.7 

2.5 

2.2 

Cross- track  roughness  angle  (8  ) 

3.8  3.5 

3.4 

3.1 

2.7 

2.5 

2.2 

Substituting  the  deep  water 
(13)  gives 

X  (8w  V2/g)  tan2  8  , 

c  c 

gravity  wave 

i  dispersion 

relation 

in  equation 

(14) 

so  that  the  change  in  apparent  wavelengths 
measurement  errors  can  be  written  as 

with 

wake 

angle  and  ship 

speed 

u=*  '=]  *vt  [ 

161V2tan  ec 

.  g 

2 

sec  8 

c 

]-. 

a 

(15) 

Using  approximate  JOWIP  parameters  (V  -  8  m/s,  8 

-  3* 

,*-o* 

) ,  we 

obtain 

41^%  0.113  AV  +  17.251  A8,  which  means  that  the  apparent  surface  wavelength 
will  be  in  error  by  approximately  30  cm  per  degree  wake  measurement  error. 
The  variation  with  ship  speed  is  much  less  significant  at  approximately  5  cm 
per  knot  where  the  average  ship  speed  error  is  estimated  at  +  0.5  knot. 

Converting  the  measured  inner  wake  angles  for  azimuth- traveling  ships 
(table  8)  into  estimated  surface  wave  wavelengths  (table  15)  by  equation  (14) 
indicates  gravity  wavelengths  considerably  in  excess  of  first  order  Bragg 
resonance,  barring  as  yet  undiscovered  bias  in  the  SAR  angle  measurements. 

As  digitally  processed  ground-range-corrected  images  like  the  example  for 
run  8-2  (figure  27)  (Shuchman,  1984)  became  available,  refined  wake  angle 
measurements  will  help  provide  better  SAR  estimates  of  surface  wavelengths  for 
comparison  with  photographic  results.  The  measured  wake  angle  from  this  image 
is  almost  1  degree  narrower  than  the  corrected  estimate  based  on  measurements 
from  the  slant  range  image . 

During  JOWXP,  we  emphasised  the  horizontally  polarized  L-band  SAR  mode. 
Use  of  the  X-band  JOWXP  images  was  severely  restricted  by  antenna  problems 
which  made  most  of  thase  images  almost  unusable  at  the  near  nadir  (23*) 
SKASAT-lika  incidence,  where  our  primary  wake  coverage  was  grouped.  Xt  must 
be  pointed  out  that  the  X-band  ship  waka  data  have  been  reported  elsewhere 
(Hammond  at  al.,  1983)  and  that  appropriate  explanations  of  wake  wave  imaging 


Figure  26.  Definition  of  radar  and  modeled  wake  angles 


For  ships  with  significant  beam,  such  as  QUAPAW  (length  60  m,  beam  12  m, 
draft  4  m) ,  the  source  of  the  shorter  wavelengths  will  be  assumed  to  be  at  the 
center  of  the  ship  so  that  a  revised  wake  angle  0^  -  a  tan  (Yo-6)/Xq)  will  be 
used  for  measurement  purposes.  This  will  reduce  the  half  wake  angles  of 
interest  in  our  nominal  1  -km  aft  wake  area  by  ~  0.3  degree  and  will  result  in 
the  AX  values  shorn  in  table  1 3  for  the  range  of  QUAPAW  wake  angles  of 
interest.  This  later  (0^)  is  the  angle  listed  for  all  photo  wake  data. 

Table  13.  SAR  versus  estimated  model  wake  angles 
for  1 -km  aft  wake  measurements. 


9r 

(deg) 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

ew 

(deg) 

1.2 

1.7 

2.2 

2.7 

3.2 

3.7 

AX 

(m) 

299 

172 

156 

127 

107 

93 

Recognizing  that  the  enhanced  populations  of  ship-generated  short  wave¬ 
lengths  identified  in  the  photographs  along  the  observed  wake  arms  propagate 
with  group  velocity  Cg,  roughly  perpendicular  to  the  ship  track,  we  may  write 

the  cross- track  roughness  radar  wake  angle  (0  )  as 

c 


0 

c 


{f) . 


(13) 


which  for  0  small  and  *  as  o  is  approximately  equal  to  the  more  general 
asin[(Cg/V)  cos  (#  -  e)]  suggested  by  Munk  (1984). 

A  comparison  between  equation  (13)  and  the  Kelvin  wake  angle  for  first  and 
second  order  Bragg  resonant  wavelengths  for  various  surface  radar  incidence 
angles  is  shown  in  table  14.  The  maximum  variation  of  Kelvin  wake  components 
from  the  ship  cross-track  direction  is  only  7  degrees  and  the  wake  angles 
differ  only  by  0.1  degree,  which  is  clearly  beyond  the  precision  of  the 
present  measurements.  Actually  we  estimate  the  maximum  error  in  the  measure¬ 
ment  of  surface  radar  wake  angles  presented  in  this  report  to  be  approximately 
±  1  degree  which  means  that  additional  work  at  higher  resolution  and  as  far 
aft  as  possible  will  be  required  to  differentiate  between  the  first  and  second 
Bragg  resonant  angles. 
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Figures  24.  CFAV  ENDEAVOUR  and  USS  QUAPAW  relative  positions  in  Dabob  Bay, 
run  8-2  (courtesy  of  DREP). 


E/O  POLAR  ANGLE  (DEGREES  CLOCKWISE  FROM  WAKE  AXIS) 


SLOPE  COMPONENT  PERPENDICULAR  TO  RADAR  BEAM 


E/O  POLAR  ANGLE  (DEGREES  CLOCKWISE  FROM  WAKE  AXIS) 

Figure  25.  Wave  slope  components  perpendicular  and  parallel  to  the  radar 
beam  (courtesy  of  DREP). 


Taking  raw  height  data  from  the  different  contractors  for  the  40  compari¬ 
son  grid  points  (shown  in  figure  22C)  allows  comparison  after  the  data  sets 
are  leveled.  Since  our  desired  output  is  relative  height  over  wavelengths 
that  are  short  compared  to  the  image  size,  the  initial  heights  were  measured 
relative  to  an  estimated  reference  plane  which  displays  an  overall  tilt  caused 
by  aircraft  attitude  and  long  wave  slope.  The  data  for  model  9's  reference 
grid  points  have  been  leveled  and  the  remaining  mean  differences  reflect  a 
vertical  bias.  The  mean  differences  between  contractors  and  between  operators 
shown  in  figure  23B  do  not  affect  relative  heights. 

Preliminary  study  of  specific  grid  points  shows  a  single  operator  repeat- 
able  height  resolution  of  4.4  _+  3.5  cm  (figure  22A)  for  this  cloudy  day.  This 
primary  factor  is  dependent  on  surface  lighting  conditions  and  altitude,  and 
must  be  statistically  evaluated  for  each  stereo  model.  No  attempt  has  been 
made  or  is  planned  to  bridge  the  surfaces  between  stereomodels.  It  is  suffi¬ 
cient  to  position  the  model  and  its  small  scale  roughness  with  respect  to 
phase  on  the  transverse  wake  waves,  wake  angle,  and  distance  aft. 

INNER  WAKE  SLOPE  OBSERVATIONS 

The  CFAV  ENDEAVOUR  carried  a  two-dimensional  laser  wave  slope  instrument 
which  provided  high  resolution,  400-Hz  data  on  changes  in  the  surface  wave 
slope  during  passage  through  a  Kelvin  wake  wave  event.  When  combined  with  the 
trisponder  positioning  data  an  accurate  determination  of  surface  roughness 
changes  within  the  wake  can  be  made  for  comparison  with  the  radar  image 
intensity. 

Figure  24  shows  the  relative  positions  of  the  ENDEAVOUR  and  QUAPAW  during 
run  8-2.  These  data  provide  the  time  ENDEAVOUR  reached  a  given  wake  angle 
(for  example  3*  @  08:22:27)  to  queue  data  from  the  slope  sensor.  Figure  25 
gives  an  example  of  DREP's  data  reduction  (Hughes,  1984)  combining  these  data 
into  slope  components  perpendicular  and  parallel  to  the  radar  plane  of  inci¬ 
dence  as  a  function  of  wake  angle.  These  data  show  maall  wave  slopes  (<  3*) 

in  the  1-3  degree  portion  of  the  wake  system  and  about  6-degree  angles  on  the 

Kelvin  cusp  waves.  Background  wave  angles  were  apparently  on  the  order  of  0.5 
degree . 

These  preliminary  data  tend  to  rule  out  specular  scattering  but  more 
detailed  analysis  is  required  before  this  can  be  confirmed. 

ANALYSIS  SUMMARY 

An  important  source  of  wake  angle  variation,  especially  when  making 
measurements  close  to  the  ship,  is  the  location  of  the  point  selected  for  the 
vertex  of  the  wake  angle.  Kelvin  wake  angles  intersect  the  ship's  track  some 
distance  ahead  of  the  ship  (Inui,  1962),  while  most  measurements  of  narrow 

azimuth,  traveling- ship  wake  angles  made  on  SAR  images  use  the  ship  itself  as 

the  vertex.  Thus  SAR  observed  wake  angles  defined  in  terms  of  measured  (Xq, 
Y  )  ,  as  shown  in  figure  26  (9  -  a  tan  [Y  /X  ]),  result  in  the  SAR  angle  (9R) 

being  larger  than  the  theoret&al  wake  angle°(9  ). 
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(A)  Comparison  between  contractors 
stereomodel  setup  parameters 
(outer  camera  orientation) 


AIRCRAFT  RELATIVE  ORIENTATION3 


CONTRACTOR 

At  (ft)’ 

AZ  (ft) 

TILT2  Au4 

it5 

A 

34.996 

-0.53 

0.87*  0.244* 

4.639* 

-1.471  * 

B 

35.00 

0.00 

0.0“  -0.017 

4.4421* 

-1 .348* 

C 

34.3 

-1.31 

2.14“  0.153“ 

4.584* 

-1.314“ 
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BY  tan  (Az/Ay)  - 

IN  MODEL 

SETUP. 

3.  FOR  CONVENTIONAL  IN-LINE  STEREO  WORK,  THE  FULL  ANGLES 
MEASURED  ARE: 


»  -  PITCH  ♦  -  ROLL  *  -  FAN  ANGLE. 

4.  THE  AM  VALUE  IS  THE  LEAD/TRAIL  BETWEEN  THE  CAMERAS 
ALONG  THE  Y-AXIS. 

5.  THE  A#  VALUE  IS  THE  CONVERGENCE  ANGLE  BETWEEN  CAMERAS 
(TOE-IN) . 

•6.  THE  Ar  VALUE  IS  THE  ROTATION  ANGLE  ABOUT  THE  VERTICAL 
Z-AXIS. 


( B)  Comparison  of  separate  measure¬ 
ments  (leveled  data  points)  for 
left  and  right  camera  reference 
grids 


COMPARISON 

CAMERA 

MEAN 

DIFFERENCE 

(cm) 
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DEVIATION 

(cm) 

BETWEEN 

LEFT 

7.0 

5.1 

CONTRACTORS 

(  A-B) 

RIGHT 

7.7 

7.4 

BETWEEN 

LEFT 

6.8 

6.2 

OPERATORS 

(C) 

RIGHT 

NO  COMPARISON 

AVAILABLE. 

(C)  Three-Dimensional  Plot: 

taken  from  spot  reading  grid 
of  2562  points,  showing  image 
position  in  wake.  (Model  9) 


(D)  Wake  profiles  near  QUAPAW 
(Model  12) 
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Figure  23.  Stereomodel  analysis  approaches. 
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Figure  22.  Stereo  spot  height  snslysis  of  wmke. 


directly  aft  of  QUAPAW  (model  12)  end  canted  at  15  degrees  to  match  the  ship 
heading  across  the  image. 


Table  12.  Stereomodels  analyzed. 


Naka  Inage  Pass  Nodal  • 

Contractor1 

Approach 

Data  Points3 

Ha  trie  boor  DBZP  7-1  1 

A 

Stereo  orientation  for  acala  (3) 

42 

B 

Not  reported 

0 

C 

Contour  plot 

2500 

Stbd  Nak*  12.6*  DUP  7-1  9 

A 

Spot  readings  on  grid 

2562 

B 

Profiles 

5473 

C 

Stereonodel  setup  only 

15 

Canter-line  £  HUP  7-1  12 

A 

Contour  plot 

50B6 

Starn  QOAPAW 

A 

Stereo  orientation  for  acala  (3) 

75 

A 

Spot  readings  on  QOAPAW  stern 

ISO 

A 

Wake  profiles  behind  QOAPAW 

260 

B 

Profiles  an  grid 

6596 

C 

Contour  plot 

40C 

C 

Spot  readings  on  grid 

300 

Port  naar  a  hip  CUP  8-3  163 

A 

Spot  readings  on  grid 

1639 

a 

Profiles 

5473 

c 

Orientations  for  study  (2) 

30 

total  points 

340004 
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grid  points  read  par  nodal, 
to  cloud  reflections, 
is:  (A)  214,  (B>  200,  and  (C)  194. 

In  order  to  heve  a  common 

data  set 

for  each  model  from  three 

contractors 

the  fixed  Reseau  grid  pattern  appearing  in  each  separate  image  was  chosen  to 
provide  common  points  of  measurement  under  which  to  determine  water  surface 
elevation.  Figure  22B  shows  the  grid  patterns  of  both  left  and  right  cameras 
as  they  appear  in  the  overlap  model.  These  reference  grid  data  from  the  three 
contractors  for  model  9  are  shown  in  figure  22C.  To  obtain  this  data  each 
company  arrived  at  their  own  set  of  spatial  parameters.  Figure  23A  gives  a 
comparison  of  three  values:  the  baseline  (Ay),  tilt  of  the  aircraft  (Az) ,  and 
left-right  orientation  angle  differences  (Aw,  A+ ,  Aie). 

The  convergence  (A+)  between  the  two  cameras  was  set  and  measured  on  the 
ground  at  5  +  0.5  degree.  The  rotation  (Ak)  was  initially  estimated  from  the 
photos  at  1.5  degrees.  The  third  angle  (Au)  measures  how  far  the  axis  of  one 
camera  leads  the  other  along  the  flight  direction.  Extraction  of  this  setting 
from  the  image  is  partially  hidden  in  forward  blur  caused  by  aircraft  motion 
during  the  minimus  exposure  (1/500th  second  ~  10  cm  on  the  ground).  Some 
variation  between  contractors  exists  for  this  set-up  value  but  this  parameter 
is  not  in  the  parallax  direction  and  does  not  affect  the  metric  solution  for 
rigidly  mounted  cameras  with  synchronized  shutters. 


These  frames  were  selected  for  analysis  in  the  absence  of  positioning 
data.  Mow  that  precise  positions  are  available  (Hughes,  1984),  this  method, 
complemented  and  compared  with  stereo  measurements,  will  be  employed  further 
back  in  the  wake  where  photographic  roughness  within  individual  SAR  resolution 
cells  can  be  compared  with  SAR  image  intensity. 

STEREO  MENSURATIONS 

Using  an  airborne  stereo  photogramme trie  approach  is  the  only  known  way  to 
extract  instantaneous  sea  surface  topography  with  fine  height  resolution  over 
even  a  moderately  large  area,  as  shown  in  Pierson  (1976).  The  need  for  pre¬ 
cise  mapping  of  small  waves  in  a  ship  wake  with  centimeter  height  resolution 
causes  a  drastic  revision  of  any  of  the  previous  mapping  approaches,  such  as 
used  by  the  Stereo  Nave  Observation  Project  (SWOP)  conducted  from  1954  to 
1960,  (Cote,  1960).  Instead  of  two  widely  separated  photo  aircraft  operating 
at  1000-2000  meters  altitude,  for  the  J0WIP  experiment,  a  unique  configuration 
was  created  for  use  at  heights  between  15-75  meters.  Rigidly  mounting  two 
identical  motor-drive  cameras  in  the  wing-tip  tank  fairings  of  a  single  air¬ 
craft,  a  fixed-baseline  system  was  created  with  a  separation  of  10.7  meters. 
With  accurate  knowledge  of  this  fixed-baseline,  we  avoided  almost  all  work  and 
problems  associated  with  pitch,  roll,  and  yaw  determination  and  ground  control 
surveying . 

As  substantiated  by  an  expert  in  fixed- baseline  photogrammetry  (Kenefick 
1984),  for  adequate  reconstruction  of  the  the  scene  in  an  analytical  stereo 
plotter  and  extraction  of  relative  height  data,  only  the  camera  separation  and 
approximate  altitude  need  be  known  in  the  rigid  stereo  camera  system  created 
by  NOSC.  The  only  ground  support  needed  is  a  calibrated  horizontal  distance 
in  a  few  photo  frames  to  serve  as  metrics  for  horizontal  scale.  During  JOWIP, 
three  such  metrics  were  used  as  "floating"  ground  control:  (1)  a  tri-spar 

buoy  photographed  on  each  major  flight  pass  with  an  8-meter  triangular 
pattern,  (2)  grid  marks  on  the  stern  helicopter  deck  of  the  CFAV  ENDEAVOUR, 
and  (3)  QUAPAW'a  existing  aircraft  operations  pattern  at  the  stem. 

Table  12  summarizes  the  Initial  set  of  data  produced  by  three  independent 
photogrammetric  contractors  using  the  highly  versatile  features  of  computer- 
controlled  analytic  stereoplotters.  Note  that  the  model  frame  numbers  are  on 
different  runs  than  previous  calm-water  examples.  Since  there  is  a  variety  of 
possible  ways  to  select,  analyze,  and  display  height  topography  from  stereo 
image  pairs  (called  stereomodels) ,  three  basic  approaches  were  used:  contour 
plot,  spot— height  reading,  and  the  directed  profile  method.  A  continuous 
analysis  method,  along  paths  of  apparent  equal  height,  produces  the  contour 
plot,  and  the  output  resembles  familiar  topographic  maps.  This  approach  is 
good  for  preliminary  study,  but  any  attempt  to  manipulate  or  level  the  data 
set  for  comparison  purposes  would  require  vast  amounts  of  computer  time  and 
multiple  assumptions.  For  this  reason  the  other  two  methods  are  more  directly 
suitable  for  analysis.  The  spot— height  method  takes  height  data  off  stereo¬ 
models  only  at  predetermined  points,  usually  on  equal  grid  spacing.  Such  an 
approach  is  shown  applied  to  model  9  in  figure  22A  and  a  resultant  three- 
dimensional  plot  can  be  seen  in  figure  23C.  The  profile  approach  is  the  most 
efficient  and  also  uses  spot  height  measurement  techniques  but  along  lines 
across  the  image  optimally  oriented  to  the  wavelength  outputs  desired.  An 
example  profile  set  is  shown  in  figure  23D,  where  cross-sectional  slices  were 


The  shift  from  shorter  wavelengths  at  small  wake  angles  to  longer  wave¬ 
lengths  at  wider  angles  is  apparent  but  complicated  by  the  additional  vari¬ 
ables  of  distance  aft  and  phase  on  the  transverse  wave.  While  the  wavelength 
bin  centered  on  30  cm  shows  maximun  enhancement  near  the  predicted  2.5-degree 
wake  angle,  the  70-cm  wavelength  bin  displays  much  higher  enhancement  for  all 
angles  from  the  wake  centerline  to  3.0  degrees.  For  larger  wake  angles,  image 
power  is  generally  depleted  in  both  these  shorter  wavelength  lines  with  en¬ 
hancement  at  the  longer  wavelengths. 

Image  power  propagation  direction  variation  shows  a  large  increase  in  the 
cross- track  direction  (90*).  Power  in  the  cross- track  bin  was  summed  over  a 
30-degree  sector  between  75  and  105  degrees  relative  to  the  QUAPAW  heading. 
Unfortunately  the  sun  azimuth  was  107  degrees  relative  at  the  time  of  this 
photo.  But  the  critical  angle  from  table  *0  shows  that  a  40-degree  slope 
would  be  required  for  direct  sun  reflection.  While  this  increase  is  not 
present  in  the  background  areas,  the  tilt  of  the  longer  transverse  waves  is 
thought  to  combine  with  an  increase  in  the  otherwise  uniform  sky  radiance  near 
the  sun,  resulting  in  a  nonlinear  enhancement.  Nonetheless,  comparison  of  the 
90-degree  bin  with  the  120-degree  bin  which  includes  the  sun  azimuth,  clearly 
indicates  a  large  increase  in  image  power  in  the  cross- track  (90*)  direction. 

In  figure  21  the  FFT  image  power  in  the  cross- track  direction  (^  15*)  is 
shown  for  the  two  shorter  wavelength  lines.  The  70  cm  plot  shows  power  spread 
along  the  3.6-degree  wake  angle  measured  on  the  SAR  image  (table  9),  but  the 
30-cm  enhancement  is  limited  to  the  area  between  200  and  300  meters  aft. 
While  these  data  clearly  indicate  large  increases  in  the  short  wavelengths 
present  along  the  narrow  wake  angles,  they  do  not  allow  us  to  select  between 
first  (30  cm)  and  second  (60  cm)  Bragg  resonant  scattering,  nor  for  that 
matter,  specular  reflection  from  the  ship  disturbed  water  surface. 
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Figure  21.  Standard  deviations  above  background  for  cross-track 
propagating  waves. 


In  this  representative  example  we  see  that  the  Image  power  has  been 
drastically  shifted  from  the  longer  to  the  shorter  wavelength  bins.  The 
maximum  image  power  now  exists  in  the  0.5-  to  0.9-meter  annulus  rather  than 
the  1.4-  to  10-meter  annulus.  The  direction  of  this  strong  trend  is  uncom¬ 
fortably  close  (120*  relative)  to  the  azimuth  of  the  sun  (107*  relative)  and 
aft  of  any  ship- induced  stationary  Kelvin  pattern  components. 


Table  11  shows  the  variation  of  total  power  (expressed  in  standard  devia¬ 
tions  above  background)  with  wake  angle  as  a  function  of  wavelength  (all 
directions)  and  propagation  direction  (all  wavelengths).  The  average  theore¬ 
tical  Kelvin  wake  angles  corresponding  to  sach  wavelength  bin  are  2.5,  3.8, 
4.9,  and  10.7  degrees,  respectively.  Kelvin  wake  wavelengths  for  the  first 
three  propagation  direction  bins  would  be  40,  30,  and  10  meters,  respectively; 
hence  showing  very  low  power  in  a  10  x  10  meter  image.  The  wavelength 
approaches  zero  in  the  cross- track  direction  where  most  of  the  power  is  shown, 
and  there  should  be  no  components  of  the  stationary  Kelvin  pattern  abaft  the 
beam. 


•  -V 


Table  11.  Total  image  power  variation  as  a  function  of  wake  angl&,  wavelength, 
and  propagation  direction  (Run  8-3) t  standard  deviations  above  background. 
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Iaage  power  in  four  directional  wave  msnber  bine  in  each  of  six  directions 
relative  to  the  QUAPAW  track  was  then  summed  and  expressed  as  a  percent  of  the 
total  image  power.  An  example  of  this  for  an  FFT  of  an  area  at  a  2.2-degree 
wake  angle,  248  meters  aft  of  the  ship,  number  115  from  photo  frame  13  is 
shown  to  the  left  in  figure  20. 
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Figure  20.  Image  power  (%),  background  (»),  and  image  (std.  dev.  above 
background),  binned  by  wavelength  (m)  and  propagation 
direction  ( deg ) . 


An  identical  procedure  was  performed  on  the  four  background  areas  shown 
and  the  resulting  statistics  (mean,  standand  deviation  on  number  of  points) 
are  shown  for  each  bin  in  the  center  of  figure  20.  Total  percents  are  also 
shorn  for  each  direction  (sector)  and  wavelength  bin  (annular  ring).  The 
background  statistics  reveal  predominant  image  power  in  the  general  direction 
(30*  relative)  of  the  surface  wind  (2.6  m/s).  The  highest  power  is  found  in 
the  1.4-  to  10-meter  wavelength  bins  as  would  be  predicted  (6.5  m)  for  this 
wind  speed. 


After  the  variations  in  total  image  power  between  FFTs  had  been  removed  by 
the  use  of  percentages,  the  percent  in  each  wavelength-direction  bin  was 
expressed  in  standard  deviations  above  or  below  background  for  comparison 
purposes.  The  diagram  on  the  right  of  figure  20  is  an  example  of  this 
procedure . 
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Figure  19.  Location  of  FFT  areas  in  wake  with  photo  example. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions  which  can  be  drawn  from  the  above  data  vary  in  confidence 
level  and  will  be  presented  here  highest  first,  followed  by  a  few  suggestions 
on  the  application  and  verification  of  these  concepts  in  future  experiments. 

The  hypothesized  ship-induced  surface  roughness  has  been  shown  to  be 
responsible  for  the  narrow  wake  angles  observed  within  the  classical  Kelvin 
wake  envelope  on  SEASAT-like  SAR  images.  These  angles  are  observed  during 
periods  of  locally  depleted  wind  wave  population:  i.e.,  calm  winds  in  the 
immediate  area  of  the  wake  pattern.  As  the  average  wind  speed  increases  above 
2.5  m/s,  narrow  wake  arm  enhancement  gives  way  to  a  full  pattern  of  transverse 
waves  against  the  wind-roughened  background  within  the  Kelvin  vee-envelope  on 
JOWIP  L-band  images.  Above  3.5  m/s  only  a  trace  of  the  turbulent  wake  is 
evident  and  above  9  m/s  all  wake  remnants  are  obscured. 

Measurement  of  ship-induced  surface  roughness  wavelengths  in  the  areas  of 
maximum  SAR  modulation  along  the  narrow  wake  arms  in  JOWIP  images  during 
periods  of  calm  winds  are  presently  insufficient  to  differentiate  between 
first  and  second  order  Bragg  scattering.  While  the  data  presented  here  show 
no  indication  of  specular  scattering,  better  statistics  are  required  to 
dismiss  this  mechanism  completely. 

Photographic  observations  on  calm  water  reveal  substantially  increased 
populations  of  short  gravity  wavelengths  along  the  narrow  angles  of  the  ship 
wake  corresponding  to  stronq  SAR  modulations  (14  dB) .  These  modulations  are 
roughly  spaced  along  track  at  the  maximum  transverse  wake  component  wavelength 
and  produce  the  largest  SAR  image  modulation  at  azimuth  ship  headings.  This 
intense  modulation  for  azimuth- traveling  transverse  waves  strongly  suggests 
that  an  orbital  motion  mechanism  like  velocity  bunching,  as  described  by 
Alpers  et  al.,  (1979,  1981,  1983),  is  the  response .  ie  mechanism. 

Further  analysis  of  transverse  wake,  wave  amplitude,  port- to- star board 
wake  arm  azimuth  offset  and  smear  may  provide  important  additional  data  for 
evaluating  SAR  image  modulation  theories.  SEASAT  observations  suggest  that 
long  ocean  swell  may  also  be  made  detectable  on  calm  seas  through  ship- 
qenerated  roughness  over  a  range  of  wavelengths  and  propagation  directions, 
but  further  analysis  of  narrow  wakes  in  locally  calm  wind  areas  on  SEASAT 
images  is  required  to  substantiate  this  hypothesis. 

The  combination  of  the  long  transverse  wave  train  and  cross-track  short 
qravity  wave  Braqo  rouqhness  that  a  ship  produces  during  calm  wind  conditions, 
provides  the  ideal  controlled  surface  conditions  essential  for  SAR  image 
modulation  experiments.  It  is  therefore  recommended  that  further  experiments 
exploiting  this  method  be  performed  to  provide  data  for  theoretical  SAR 
modelinq.  In  each  such  experiment,  after  the  baseline  case  (displaying  a 
narrow  transverse  wave  modulated  wake  for  an  azimuth  traveling  ship  on  calm 
water)  has  been  demonstrated,  the  effects  of  transverse  wave  direction  (ship 
heading),  wavelenoth  (ship  speed),  and  incidence/scattering  mechanism  (ship 
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locations  in  the  swath)  under  increasing  background  wind  conditions  can  be 
quantitatively  evaluated.  Images  collected  in  this  manner,  when  combined  with 
simultaneous  photographic  records  of  surface  conditions  and  ambient  wind  and 
surface  wave  spectral  measurements,  can  provide  the  data  essential  for 
quantitative  comparison  of  SAR  ocean  wave  imaging  theories. 


REFERENCES 


Alpers,  w.R.  and  C.L.  Rufenach.  "The  Effects  of  Orbital  Motions  on  Synthetic 
Aperture  Radar  imaging  of  Ocean  Waves.”  IEEE  Transactions  on  Antennas 
Propagation,  1979,  AP-27,  685-690. 

Alpers,  W.R. ,  D.B.  Ross,  and  C.L.  Rufenach.  ”Cn  the  Detectability  of  Ocean 
Surface  Waves  by  Real  and  Synthetic  Aperture  Radar,”  journal  of  Geophysical 
Research,  1981,  86,  6481-6498. 

Alpers,  w.R.  "Monte  Carlo  Simulations  for  Studying  the  Relationship  Between 
Ocean  Wave  and  Synthetic  Aperture  image  Spectra.”  journal  Geophysical 
Research,  1983,  88,  1745-1759. 

Cote,  L.  J. ,  J.O.  Davis,  W.  Marks,  R.J.  McGough,  E.  Mehr,  W.J.  Pierson,  J.F. 
Kopek,  G.  Stephenson  and  R.C.  Vetter.  "The  Directional  Spectrum  of  a  Wind 
Generated  Sea  as  Determined  from  Data  Obtained  by  the  stereo  Wave  Observation 
Project."  Meteorol.  Paper,  Vol.  2,  No.  6,  New  York  University  press.  New  York, 
June  1960. 

Cox,  C.,  and  W.H.  Munk.  "Statistics  of  the  Sea  Surface  Derived  from  Sun 
Glitter,"  journal  Marine  Research,  1954,  13(2),  198-227. 

Fu,  L.L. ,  and  B.  Holt.  " SEASAT  views  Oceans  and  Sea  Ice  with  Synthetic 
Aperture  Radar,"  Publication  81-120,  Jet  Propulsion  Laboratory,  Pasadena,  CA, 
Feb  15,  1982. 

Hager,  R.o.  "The  side-looking  radar  image  of  time-variant  scenes,"  Radio 
Science,  1980,  J_5,  (4),  749-756. 

Havelock,  T.H.  "Wave  Patterns  and  Wave  Resistance,"  irans.  INA,  1934  74, 
430-442. 

Hughes,  p.A.  private  Communication,  "internal  Waves  in  Dabob  Bay  During 
JCWIP",  7  Mov  1983. 

Hughes,  B. A.  Private  Communication,  " QUA PAW  Speeds,  Headings  and  ENDEAVOUR 
Wake  Crossing  (3°)  Times",  25  Apr  1984. 

Hulbert,  F.O.  "The  Polarization  of  Light  at  sea,"  Journal  of  the  Optical 
Society  of  America,  1934,  2A,  (35). 

Inui ,  T.  "Wave-Making  Resistance  of  Ships,"  Trans.  Soc.  Nav.  Arch.  Ehgr., 
1962,  70,  283. 


67 


Keller,  W.C.  and  J.W.  Wright.  “Microwave  Scattering  and  the  Straining  of 
Wind-Generated  Waves Radio  Sci.,  1975,  _1_0*  139-1947. 

Kenefick,  J. P.  "Analytical  Study  of  Ocean  Waves  and  Ship  Wakes  Using  Airborne 
Stereo  Photogramme try."  Letter  Report,  Contract  No.  N66001 -84-M-01 1 3,  6  June 
1984. 

Lighthill,  J.  "Waves  in  Fluids,"  New  York,  Cambridge  University  Press,  1978. 

Munk,  W. H.  "Kelvin-Wake-Dabob  Bay,"  in  JASON  SEASAT  III  and  IV,  Aug  1984, 
Mitre  Cbrp.  JSR-84-203,  3-1  to  3-15. 

Pierson,  w.  J.  "The  Theory  and  Applications  of  Ocean  Waves  Measuring  Systems 
at  and  Below  the  Sea  Surface  on  the  Land,  from  Aircraft  and  from  Spacecraft," 
NASA  Itept.  CR-26-2646,  Contract  No.  NAS  5-20041,  Jan  1976. 

Phillips,  OM.M.  "The  Dynamics  of  the  Upper  Ocean,"  51-54,  1977,  New  York. 

Cabridge  Uhiversity  Press. 

Shuchman,  R.A.,  and  J.s.  Zelenka.  "processing  of  Ocean  Wave  Data  from 
Synthetic  Aperture  Radar,"  Boundary  Layer  Meteorol,  1978,  1 3,  181-191. 

Shuchman,  R.A.,  "Processing  of  Synthetic  Aperture  Radar  Data  of  Ocean  Waves, 
in  Oceanography  from  Space."  Pd.  J. F.R.  Gower,  477-496,  1981,  New  York, 
Plenun. 

Shuchman,  R.A.  "Digital  Ground  Range  Corrected  SAR  image,"  JOWIP  Run  8-2,  June 
1984. 

Stillwell,  Jr.,  D.  "Directional  Fnergy  spectra  of  the  Sea  from  Photographs ," 
Journal  of  Geophysical  Research,  1969,  J74,  1974-1986. 

Valenzuela,  G. R.  "Theories  for  the  Interaction  of  Flectromagnetic  and  Ocean 
Waves  -  A  Review,"  Boundary  layer  Meteorol.,  1978,  _1_3,  61-85. 

Vesecky,  J. F. ,  and  R.H.  Stewart.  "The  Observation  of  Ocean  Surface  Phenomena 
Using  Imagery  from  the  SEASAT  Synthetic  Aperture  Radar,"  J.  Qeophys.  Res., 
1982,  87,  3397-3430. 

Vesecky,  J. F. ,  and  R.H.  Stewart,  R.  A.  Shuchman,  H.M.  Assas,  F.R.  Kasischke, 
and  J.D.  Lyden.  "Ch  the  Ability  of  Synthetic  Aperture  Radar  to  Measure  ocean 
Waves,"  Proceedings  of  the  IUCRM  Symposium  on  Wave  Dynamics  and  Radio  Probing 
of  the  Ocean  Surface,  O.M.  Phillips,  editor.  Sept  1982. 

Webster,  W.D.  "Private  Discussions  on  Ship  Wake  Structures,"  Dept,  of  Naval 
Architecture,  Uhiv.  of  Calif.,  Berkeley,  5  Aug  1982. 


Wright,  J.W.  "A  New  Model  for  Radar  Sea  Clutter,”  ZBEE  Transactions  on 
Antennas  and  Propagation,  1968  AP-1 6,  217-223* 

Wright,  J.W.,  w.j.  Plant,  and  W.C.  Keller.  "Ocean  Wave-Radar  Modulation 
transfer  Functions  from  the  West  Coast  Bcperiment,”  journal  of  Geophysical 
Research,  1980,  85  4957-4966* 


APPENDIX  A  -  AVAILABLE  DIGITAL  DATA 


There  are  three  categories  of  digital  data;  tiae  series, 
image  —  all  are  recorded  on  9- track  1600  BPI  digital  tape. 


track ,  and 


The  time  aeries  data  and  track  are  both  stored  in  VAX/VMS  FILES-II  level 
2  direct  access  files  and  nay  be  read  from  tape  using  the  VMS  "COPY-  command. 
The  Images  are  stored  integers  with  one  row  per  record. 
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A-2  TIME  SERIES  DATA 


mwTO 

a.  Matar  valocity  Measured  with  raspact  to  tha  ship  expressed  as  North-  and 
East-going  coaponants 

b.  Distanca  tha  ship  has  traveled  through  the  water;  obtained  by  integrating 
the  bow  current  aeter  signal  (asters) 

c.  Conductivity,  teaperatu re,  sigaa-t,  and  depth  below  the  mean  sea  surface 
(Siaasns/aeter)  (Degrees  C)  (Kg /a3)  (asters) 

d.  Tine-differentiated  wave  height  (awtera/second) 

e.  Ships  heading  (degrees) 

f.  Relative  wave  height  (asters) 

g.  Mean-square  differentiated  wave  height  ((m/s)3) 

h.  Mean-square  wave  height  (aeters  squared) 

i.  Mean-square  wave-slope 

)•  Sane  as  a.  but  aaasured  free;  stern  current  aster 

k.  Sane  as  b.  but  aaasured  froa  stern  current  aeter 

l.  Wave-slope,  north  and  east  coaponants  with  laser  intensity  for 
each  pair  of  coaponants 

a.  Wind  valocity  (a/s) 

Tapes  are  available  for  the  following  SAR  runs. 

5-3,  5-4,  6-4,  8-1,  8-2,  8-3,  8-4,  9-1,  9-2,  and  9-3 


A- 3  SAR  IMAGES 
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APPENDIX 

»  -  VERTICAL  (STEREO) 

PHOTO  TIMES 
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lime  (PDT) 

10:  16:29 

1 

11 :01 :04 

1 

15: 18:  36 

10:16:30 

2 

11:01:05 

2 

1 5:18:37 

10:  16:31 

3 

1 1 :  01  :  06 

15:18:38 

10:16:32 

4 

1 1:01 :07 

4 

15:18:39 

10: 16:33 

5 

1 1  :  01  :  08 

5 

15: 18:40 

10:16:34 

6 

1 1:01 :09 

6 

15:18:41 

10: 16:35 

7 

1 1  :  01  : 10 

7 

1 5: 18:42 

10:16:36 

8 

11:01:11 

8 

15:18:43 

1  0: 1 6:  37 

9 

11:01:12 

9 

15:  18:44 

10:16:38 

10 

11:01:13 

10 

15:18:45 

11 

1 1 :  01 : 1 4 

11 

15: 18:46 

12 

11:01:15 

12 

15:18:47 

13 

11:01:16 

13 

1  5: 18:48 

14 

1 1:01:17 

14 

15:18:49 

15 

11:01:18 

15 

15: 18:50 

16 

11:01:19 

16 

15:18:51 

17 

15:18:52 

Run  2 

18 

15: 18:53 

19 

15:18:54 

Fr# 

lime  (PDT) 

20 

15:18:55 

1 

11 :05: 15 

2 

1  1  :05:16 

3 

11 :05: 17 

4 

1  1  :05:18 

5 

11:05:19 

6 

1 1:05:20 

7 

11:05: 21 

8 

1 1:05:22 

9 

11:05: 23 

10 

1 1 :05: 24 

11 

11:05: 25 

VJS* 


.VvIn^vvVIv'Is^v'jw’Iv' 


**  .*•  ;**  ^  ;•►>*■ 


DREP  6-4 
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9:46: 32 

9:46:33 

9:46:34 

9:46:35 

9:46:36 

9:46:37 

9:46:38 

9:46:39 

9:46:40 

9:46:41 

9:46:42 

9:46:43 

9:  46: 44 

9:46:45 

9:46:46 

9:  46: 47 

9:46:48 

9:46:49 

9:46:50 

9:46:51 

9:46:52 

9:46:53 

9:  46: 54 

9:46:55 

9:  46: 56 

9:46:57 

9:46:58 

9:46:59 

9:47:00 
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ROLL  NUMBER  8 
RUN  7-1 


RUN  7-3 


ROLL  NUMBER  9 


RUN  7-5 


RUN  7-7 


16:08:56 
15; 12:08 
16:12:51 
16: 15:20 
16:19:37 
17:11:92 
17:12:57 
17: 13:^0 
17: 15:37 
17:17:25 
17:19:39 

17:20:98 
17:22:55 
17:25:22 
18: 10:18 
18:15:30 
18: 18:95 
18:22:57 
18:23:39 
18:29:20 
18:25:11 
18:28:35 
13:31:91 
18:33:17 
18:39:09 
18:35:57 
16:99:20 
19:91:19 
19:92:00 
19:92:53 


16:09:57 
16: 12:18 
16:12:58 
16:16:28 
16:20:13 
17: 12:05 
17:13:05 
17:19:25 
17:16:30 
17:17:36 
17: 19:99 
17:20:57 
17:23:05 
17:25:30 
1 8 : 12:9? 
18:16:08 
18:22:09 
18:23:05 
18:23:96 
18:29:35 
18:26:18 
18:29:05 
18:32:05 
18:33:23 
18:39:37 
18:37:13 
18:97:59 
19:91:29 
19:92:06 
19:93:02 


16:10:09 
16:12:26 
16:13:51 
16:16:91 
16:20:97 
17:12:23 
17:13:15 
17:19:39 
17:16:36 
17:18:28 
17:20: 10 
17:21: 17 
17:29:05 
17:25:38 
1 8: 12:55 
18:16:52 
18:22:21 
18:23:12 
18:23:59 
18:29:93 
18:27:06 
18:30:50 
18:32:39 
18:33:33 
18:39:59 
18:93:53 
18:98:09 
19:91:31 
19:92:19 


1b: 10:53 
16:12:35 
16:19:37 
16:17:39 
16:21:16 
17:12:36 
17:13:26 
17:19:52 
17:16:56 
17:18:50 
17:20:35 
17:21:35 
17:29; 19 
18:08:91 
18:13:16 
18:18:31 
18:22:37 
18:23:20 
18:29:01 
18:25:03 
18:28:00 
18:31:10 
13:32:90 
18:33:91 
18:35:10 
18:99:00 
19:39:97 
19:91:91 
19:92:29 


16:11:25 
16:12:93 
16:15:19 
16:19:08 
17:10:99 
17:12:99 
17:13:32 
17:15:00 
17:17:07 
17: 19: 15 
17:20:90 
17:21:55 
17:25:16 
18:10:09 
18:13:98 
18: 18:90 
18:22:98 
18:23:27 
18:29:10 
18:25:90 
18:23:16 
18:31:18 
18:32:99 
18:33:56 
18:35:25 
18:99:08 
19:91:09 
19:91:59 
19:92:92 


t  ••.'■•r- * 

IJULY  28l 

7:06:52 

7:09:33 

7:12:48 

7:23:26 

7:24:29 

ROLL  NUMBER  10 

7:25:36 

7:25:47 

7:26:04 

7:26:52 

7:27:00 

RUN  8-1 

7:27:17 

7:27:23 

7:27:31 

7:27:39 

7:27:48 

7:28:01 

7:28:32 

7:29:41 

7:32:41 

7:33:12 

7:34:41 

7:35:06 

7:35:23 

7:35:37 

7:40:12 

7:41:13  P 

8:04: 17» 

8: 10:07 

8:11:43 

8:11:53 

5: 12:02 

8: 12:10 

8:12:19 

8:12:28 

8:12:38 

8: 12:48 

8: 12:59 

8: 13:10 

8:13:21 

8:15:07 

8:15:33 

8:16:31 

8:17:04 

8:17:33 

8:18:48 

RUN  8-2 

8: 19:09 

8: 19:24 

8:19:33 

8:19:43 

8:19:59 

8:20:10 

8:20:32 

8:20:57 

8:21:09 

8:21:41 

8:22:33 

8:22:48 

8:22:58 

8:24:03 

8:24:59 

8:56:36 

8:56:46 

8:57:24 

8:57:51 

8:57:59 

8:58:05 

ROLL  NUMBER  11 

8:59:14 

8:59:21 

8 : 3  j ;  3 1 

8:59:39 

8:59:56 

9:00:02 

9:00:40 

9:00:49 

9:01:01 

9:01:16 

RUN  8-3 

9:02:22 

9:02:39 

9:02:48 

9:02:55 

9:03:10 

9:03:32 

9:03:41 

9:03:50 

9:04:08 

9:04:43 

9:05:06 

9:05:14 

9:05:24 

9:05:31 

9:05:41 

9:06:07 

9:06:14 

9:06:21 

9:06:37 

9:06:46 

9:06:54 

9:07:04 

9:07:12 

9:07:20 

9:07:29 

9:07:48 

9:08:19 

9:08:42  T 

9:49:08 

9:49:14 

9:49:24 

9: 50:06 

9:50:19 

9:50:27 

9:50:33 

9:50:40 

9:50:4? 

9:50:56 

9:51:08 

9:51:16 

9:51:23 

9:51:30 

9:51:36 

9:51:43 

9:51:50 

RUN  8-4 

9:51:58 

9:52:07 

9:52:14 

9:52:26 

9:52:58 

9:53:06 

9:53: 15 

9:53:21 

9:53:47 

9:53:54 

9:54:01 

9:54:08 

9:54:15 

9:54:31 

9:54:52 

9:54:57 

9:55:05 

9:55:17 

9:55:26 

9:55:35 

9:55:45 

9:55:52 

9:56:00 

9:56:16 

9:56:31 

ROLL  NUMBER  12 

13:09:33 

13:09:49 

13:10:43 

13:10:51 

13:11:21 

13:11:40 

13:12:47 

13:12:55 

13:13:48 

13:13:59 

13:14:19 

13:14:32 

13: 14:48 

13:15:00 

13:15:22 

RUN  9-1 

13:15:35 

13:15:57 

13:16:38 

13:16:45 

13:16:57 

13:17:13 

13: 17:24 

13: 17:31 

13:17:57 

13:18:07 

13:18:18 

13:18:30 

13:18:45 

13:18:56 

13:19:49 

13:20:17 

13:20:37 

13:21:02 

13:22:02  T 

13:53:64 

13:53:15 

13:53:52 

13:54:21 

13:55:23 

13:55:32 

1  55:41 

13:55:52 

13:56:00 

13:56:06 

13:56:11 

13:56:15 

13:56:20 

13:56:27 

13:56:35 

13:56:52 

RUN  9-2 

13:57:04 

13:57:11 

13:57:22 

13:57:55 

13:58:05 

13:58:28 

13:58:39 

13:58:46 

13:58:54 

13:59:01 

13:59:09 

13:59:16 

13:59:37 

13:59:49 

13:59:56 

14:00:00 

14:00:05 

14:00:11 

14:00:40 

14:00:48 

14:00:57 

14:01:10 

14:01:19 

14:01:28 

14:01:37 

14:01:51 

14:02:14 

14:02:22 

14:02:32 

ROLL  NUMBER  13 

14:37:48 

14:38:56 

14:39:08 

14:39:15 

14:39:21 

14:39:46 

14:39:58 

14:40:12 

14:40:41 

14:40:48 

14:40:56 

14:41:07 

14:41:35 

14:41:40 

14:41:48 

14:42:13 

14:42:19 

14:42:26 

14:42:37 

14:43:18 

14:43:24 

14:43:42 

14:43:48 

14:44:05 

14:44:11 

RUN  9-3 

14:44:25 

14:44:31 

14:44:37 

14:45:27 

14:45:34 

14:46:05 

14:46: 10 

14:46:16 

14:46:21 

14:46:26 

14:48:50 

14:49:00 

14:49:05 

14:49:18 

14:46:30 

14:46:35 

14:46:40 

14:46:49 

14:46:55 

14:47:02 

14:47:09 

14:47:15 

14:47:21 

14:47:26 

14:47:40 

14:47:45 

14:47:55 

14:48:00 

14:48:07 

14:48:12 

14:48:18 

14:48:26 

14:48:33 

14:48:42 

14:49:35 

14:49:44 

14:49:51 

14:50:01 

14:50:07 

14:50:16 

14:50:27 

14:50:36 

14:51:03 

14:51:08 

14:51:13 

14:51:21 

14:51:36 

14:51:41 

14:51:47 

14:51:56 

14:52:05 

14:52:38 

14:52:44 

80 


f  T* 


